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Abstract 

We study the Cauchy problem for general, nonlinear, strictly hyperbolic systems 
of partial differential equations in one space variable. First, we re- visit the construc- 
tion of the solution to the Riemann problem and introduce the notion of a nondegen- 
erate (ND) system. This is the optimal condition guaranteeing, as we show it, that 
the Riemann problem can be solved with finitely many waves, only; we establish that 
the ND condition is generic in the sense of Baire (for the Whitney topology) , so that 
any system can be approached by a ND system. Second, we introduce the concept of 
inner speed variation and we derive new interaction estimates on wave speeds. Third, 
we design a wave front tracking scheme and establish its strong convergence to the 
entropy solution of the Cauchy problem; this provides a new existence proof as well 
as an approximation algorithm. As an application, we investigate the time-regularity 
of the graph solutions (X, U) introduced by the second author, and propose a geo- 
metric version of our scheme; in turn, the spatial component X of a graph solution 
can be chosen to be continuous in both time and space, while its component U is 
continuous in space and has bounded variation in time. 
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1 Introduction 



In this paper, we are interested in entropy solutions to general, nonlinear, strictly hyper- 
bolic systems of partial differential equations in one space variable, and we investigate 
several important issues of the theory, especially the Riemann problem, the wave interac- 
tion estimates, and the graph solutions. In particular, a version of the wave front tracking 
scheme is introduced below, which extends schemes originally introduced by Dafermos 
[TJ], DiPerna [IT], Bressan [9], and Risebro jlOj for genuinely nonlinear systems of 
conservation laws (see also [1]). The new version applies to general systems that need not 
be genuinely nonlinear nor in conservative form. 

Recall that the front tracking scheme is a variant of the random choice scheme intro- 
duced by Glimm [PS] for genuinely nonlinear, conservative systems. Recently, Glimm's 
scheme was extended to conservative systems with general flux by Bianchini [7J, Iguchi 
and LeFloch [23], and Liu and Yang [38]. The key contributions in the above works 
were the derivation of sharp estimates for the solution of the Riemann problem, together 
with the introduction of suitable interaction functionals controling the total variation of 
solutions to the Cauchy problem. The analysis in [7J E31 EE] took its roots in earlier works, 
pioneered by Liu in 1981, on particular systems (see Liu [SI EH1 EE1 E7J , LeFloch et 
al. [2SI EH [30], Chern p3], and Ancona and Marson [H El E]), and in the recent 
developments on the vanishing viscosity method by Bianchini and Bressan [8]. 

In comparaison to the Glimm scheme, implementing a front tracking scheme is more 
demanding since, in addition to the interaction estimates and functionals controling the 
growth of wave strenths (provided in [7J E31 EH] ) , one also needs a precise control of wave 
speeds. In a front tracking method, when dealing with a general hyperbolic system and in 
order to guarantee that the limiting function is the entropy solution of interest, it is often 
necessary (at interactions) to split a given front into several fronts (of smaller strength, 
say). Clearly, it is preferable to avoid to repeat this splitting step too often during the 
evolution — for otherwise this could lead to a finite-time blow-up of the algorithm with 
the number of wave fronts or of interactions points becoming infinite. 

We found it convenient to introduce wave fronts that are not just single shocks or single 
rarefactions but, rather, wave packets consisting of (finitely or possibly infinitely) many 
shock and rarefaction waves propagating at the same speed. To control the speed of a j- 
wave packet (u_, u + ) we introduce the notion of inner speed variation $,•(?/_, u + ), defined as 
the largest minus the smallest wave speeds within the wave packet under consideration. One 
contribution in the present work is to provide sharp interaction estimates on i9j(u-,u + ), 
via suitable convexity and wave interaction arguments. 

The present paper also provides a significant generalization of the earlier analysis by 
Iguchi and LeFloch [23], who dealt with general systems approached by piecewise 
genuinely nonlinear (PGNL) ones. We introduce below the notion of a nondegenerate (ND) 
system and we show that the ND condition is the optimal condition guaranteeing that the 
Riemann problem is solvable with finitely many waves, only. Moreover, the nondegeneracy 
condition is shown to be fully generic in the sense of Baire (in the Whitney topology), so 
that any system can be approached by a nondegenerate one. 
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We also include below a discussion of hyperbolic systems in nonconservative form. Re- 
call that, for such systems, the standard notion of a distributional solution does not apply, 
and it was recognized in the 90's by LeFloch (2E1 E7] and Dal Maso, LeFloch, and 
Murat [16], that one must supplement the hyperbolic system with a prescribed family 
of continuous paths, related with the "interior structure" of discontinuities in solutions 
and typically constructed from the set of all traveling wave solutions associated with a 
given parabolic regularization. The Riemann problem was solved in [TB] for nonconser- 
vative, genuinely nonlinear systems, and a suitable generalization to the Glimm scheme 
was proposed by LeFloch and Liu [33]. We observe in this paper that both techniques 
developed in [8] and [23] allow for generalizations to nonconservative systems. 

In summary, our results provide a new proof of the global existence of the entropy solu- 
tion (in the class of functions with tame variation) to the Cauchy problem associated with 
a strictly hyperbolic system, and represent an alternative approach as well as a generaliza- 
tion to the proofs given in [TBI 133] (for genuinely nonlinear systems) and, more recently, 

in PUSES EE! ■ 

On the other hand, recall that the uniqueness of the entropy solution was established 
earlier by Bressan and LeFloch [TT] (for genuinely nonlinear, conservative systems), 
Baiti, LeFloch, and Piccoli [S] (for general conservative or nonconservative systems 
and arbitrary jump relations), and Bianchini and Bressan [8] (for vanishing viscosity 
limits based the identity viscosity matrix). 

Finally, we apply our front tracking scheme and investigate the time-regularity of the 
graph solutions (X, U) introduced in LeFloch [3H [32] • We propose a geometric version 
of our scheme and, in turn, provide here a new, more regular parametrization of the graph 
solutions, which is not only continuous in space but also continuous in time, except at 
countably many times where wave cancellations/interactions take place. 

A brief outline of this paper follows. Section [2] deals principally with the Riemann 
problem and generalizes the approach developed by Iguchi and LeFloch's; in particular, 
it includes the discussion of the nondegeneracy condition, wave interaction estimates, and 
inner speed variation estimates. Section [3] provides a second proof of the inner speed vari- 
ation estimates which follows the approach developed by Bianchini and Bressan. Section H] 
contains a brief discussion of generalizations of our results to nonconservative systems. 
Next, Section [5] provides our new version of the front tracking scheme, together with the 
convergence proof. Finally, Section [U] shows how to apply the previous framework in order 
to investigate the time-regularity of graph solutions. 



2 Non-degenerate hyperbolic systems of conservation 
laws 

2.1 Notation 

In the present section and in the following one we are primarily interested in conservative 
systems. However, as we will explain in Section H] many of the forthcoming arguments 
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carry over to nonconservative systems, with only minor modifications. We thus consider a 
strictly hyperbolic system of conservation laws, 

d t u + d x f{u) =0, u = u{t,x)eR N , t>0,xeR, (2.1) 

where, as is usual, we assume that all solutions take values in a neighborhood of some 
constant state in ~R N - which is normalized to be the origin. We denote by Bs the open 
ball in M. N centered at the origin and having radius 5 > 0, and we assume that all values 
u(t,x) G Bg 1 for some 8\ > 0. In (12. ip . the map / : B$ l — > M, N is smooth and, for each 
u G Bs 1 , A{u) := Df(u) admits N real and distinct eigenvalues Xi(u) < ■ ■ ■ < Xn(u). We 
denote by rj(u) and lj(u) the left- and right-eigenvectors associated with \j(u), i.e., 

A{u) Tj{u) = Xj(u) Tj(u), lj(u) ■ A(u) = Xj(u) lj(u), 

and normalized so that 

|r,-(u)|=l, lj(u) ■ Tj{u) = 1, 
k(u) ■ Tj{u) = for i ^ j. 

In [23], the authors advocated to use, for each j-wave family, a foliation of Bg 1 based 
on a global parameter u i— > ftj(u), satisfying 

VfXj(u) ■ Tj{u) 7^ 0. 

Each wave curve s t— > ipj{ m ) = ipj{ m ] u o) is parametrized to ensure that the state ipj(m) 
lies on the submanifold {/^(it) = m|, i.e., 

^•(^•(m)) = m. 

In the present paper we adopt a particular global parametrization. We fix a constant vector 
I such that I ■ rj (u) never vanishes, 

and we set 

Hj{u) = /i(u) :— I ■ u. 

Hence, our global parameter is also independent of the specific wave family under consid- 
eration. The following formulas will involve the renormalized eigenvectors 

~ r j r j 

Recall that, for the j-integral curve 

Oj{uJ) = \ Wj(m; u_) / m+ <m< m*} 

defined by 

d m Wj = rj(vjj), Wj{m_]uJ) = 
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where m_ := /i(w_), one can check that there exists 62 < Si such that the following property 
holds. For all w_ G B$ 2 , the map [m*, m*] 3 m 1— > Wj(m_; uJ) is smooth and takes its values 
in Bg 2 , and the end points m*, m* (depending on u_) satisfy Wj{m^ w_), Wj(m*; u_) G <9£>«5 2 . 
Moreover, one has 

Wj(m; Uo) = u + (m — m ) rj(u ) + C(m — mo) 2 

with m := fi(uo). 

Similarly, the j-Hugoniot curve 

TCj(u-) := {^(m; u„) / m* < m < m*}, 

defined by 

—Xj(u-, Vj(m; u-)) (yj(m; u_) — it-) + f(vj(m] «_)) — f{uJ) = 0, 

satisfies 

Uj-(m; ito) = u + (m — m ) r^-(iio) + C( m — m o) 2 ; 
whereas the j-shock speed \j{m;uJ) := Aj(tt_, Uj(m; u_)) satisfies 

Aj(m; m ) = Aj(m ) + ^(VAj • r})(u ) (m - m ) + C(m - m ) 2 . 

As usual, we restrict attention to propagating discontinuities («_,«+) with u + G 
Hj(u-), satisfying the following entropy criterion [391 ED EU [351 ES, [37] . Setting 
m_ := /i(w_), m + := /i(w+), and 

M + = Vj(m + ; U-), 

the shock (w_,tt + ) is said to be admissible if and only if 

\j(it-,u+) < Xj(m;U-) for all m between m_ and m+. 

In other words, the shock speed achieves its minimum value at the point m+. Following 
[2H (37J [23], we will be interested in combining rarefaction curves and shock curves and 
constructing the j-wave curve issuing from a given left-hand state uq G £><j 2 , 

Wj(u ) := {iJ}j(m;uo) / m* < m < m*}, 

where the end points of the curve satisfy 

m* = mj*(u ) < h{uq) < m*(uo) = m*. 
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2.2 Riemann problem for nondegenerate systems 



Recall that, following Liu's pioneering contribution to the subject (Liu [HH [35], [36l [37]). 
Iguchi and LeFloch [23] constructed and investigated the regularity of the entropy 
solution to the Riemann problem, for systems with piecewise genuinely nonlinear (PGNL) 
flux (see below for a definition). Here, we introduce a larger class and we establish that all 
of the results in [23J remain valid. Recall that all of the key estimates depend on the C 2 
norm of the flux only, and allowed the authors in [23] to cover any flux realized as a limit 
of PGNL functions. Our more general condition will be shown to be fully generic, so that, 
by density, we can cover systems with fully arbitrary flux. 

It will be convenient to introduce the notation (1 < j < N, u G Bs x ) 



together with the following definition. Note that the definition clearly makes sense for a 
general matrix-valued mapping A which need not be a Jacobian matrix. 

Definition 2.1. The matrix-valued mapping A = A(u) (u G Bs l ) is called nondegenerate 

(ND, in short) if for every j = 1, . . . , N and every u G Bs 1 , not all of the values n\ f \u) 
(1 < k < N + 1) are zero. 



When the ND condition holds, we will also say that (12 .11) is a nondegenerate system or 
that, in the conservative case, / is a nondegenerate flux. 

For instance, in the scalar case N — 1 a function / : R — > M. is nondegenerate if and 
only if /" and do not vanish simultaneously. 

Given any nondegenerate flux / we can associate to each point u G Bs 1 its critical 
exponent p = Pj{u) as the smallest index k > 1 such that iTj k \u) ^ 0. Our definition 

includes the genuinely nonlinear flux for which ir^^u) never vanishes (so that Pj{u) = 1), 
as well as the piecewise genuinely nonlinear (PGNL) flux introduced in |23j where it 
was assumed that 7r^\u), TTj 2 \u) never vanish simultaneously (so that Pj{u) equals 1 or 2). 

The case of linearly degenerate fields for which tt^ (u) vanishes identically could be easily 
included in the present discussion, but for simplicity in the presentation we prefer to cover 
this case latter via a general density argument. 

Our first objective is to generalize the construction in [23] and to prove that the Riemann 
problem for ND flux can be solved uniquely and that the solution contains finitely many 
waves, only. 

Theorem 2.1. (Riemann problem for nondegenerate flux.) 

Let (12. ip be a strictly hyperbolic system of conservation laws with nondegenerate flux / : 
Bs 1 — > 1R JV . Then there exists 62 < Si depending only upon the C 2 norm of f , such that the 
following properties holds. 




) := rj(u) ■ VAj(m), 
(u) := Tj{u) ■ Wnf\u), 



k = l,2,..., 



(nf\u),nf\u) 



•••,< +1) H)^(o,o,...,o). 
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1. For all uq G B$ 2 and all j = 1, . . . , N there exist m* = m JVr (M ) and m* = m*(uo) 
and a mapping ipj = ipj(m; Uq) defined for m* < m < m* and Lipschitz continuous in both 
arguments, such that 

if}j(m; Uq) G Bs 2 , m± < m < m*, 

V>,-(a*(«o);«o) = u , 
^•(m*;u ), ij)j(m*;uo) G dBs 2 . 

2. For each m G [m*, m*] there exist an integer q and a finite sequence of states 
U\, u 2 , • • • , M2g+i G £><5 2 swc/i that 

and (for all relevant values of k) u<ik is connected to u 2 k+i by a j -rarefaction wave while 
U2k+i is connected to u 2 k+2 by an admissible j -shock, with 

Aj(tto) < Aj(wi) <Xj(ui, u 2 ) < Xj(u 2 ) < ■ ■ ■ 

■ • • < Xj(u 2q -i, U 2q ) < \j(u 2q ) < \j(u 2q+ l). 

3. For any Riemann data ui, u r G B$ 2 the Riemann problem (12.11) . 

/ n \ jui, x < 0, 
I « r , x > 0, 

admits a unique, self-similar solution made of finitely many rarefaction waves and admis- 
sible shock waves. 

To establish these results we will explain how the arguments in [23] can be extended 
to cover nondegenerate flux. Certain monotonicity properties of the shock speed along the 
Hugoniot curve and its critical values play a central role in the construction of the wave 
curves, and this will be the subject of the following section. 

Once the Riemann problem is solved and provided finitely many waves only arise in the 
construction, one can follow [23] and derive additional regularity and interaction estimates. 
Indeed, we emphasize that the proofs therein rely only on the facts that the Riemann 
solution contain finitely many waves and that all Hugoniot and integral curves associated 
with ( 12. ip are smooth maps. We summarize here the results that we will need later. 

Proposition 2.1. (Regularity of the wave curves.) Under the assumptions of Theorem \2.1\ 

for each j = 1, . . . , N the map ifjj = ijjj(m; u) is Lipschitz continuous with respect to both 
m,u (with Lipschitz constant depending only on the C 2 norm of f). Moreover, the first- 
order derivatives ofipj with respect to m, u are Lipschitz continuous at the point m = n(u), 
that is 

ipj(m; u) = Tj(u) + 0(m — fx(u)), 

D u ijjj(m; u) = Id — rj(u) <g> I + 0(m — /u(it)), 

where ® denotes the tensor product of two vectors. 
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To any left-hand state w_ and a right-hand state u + = i[)j(m + ;U-), we associate the 
wave strength 

6j(u-,u+) := l^{u-) — 

and the associated wave speed function Xj(s, m + ; w_), which by definition yields the 
speed of waves making up the wave fan. Given two j-waves (u-,u + ) and (u'_,u' + ) we 
introduce the generalized angle between them, as follows 



h- r m '+ (\j(s',m' + ;u'_) - \j(s,m + ;u-))_ 

Oj{u-,u + ;u'_,u' + ) : = / / — dsds', 

[m + — m_) (m + — m_) 



where m± := fi(u±) and m' ± := fi(u'±). 

Proposition 2.2. (Interaction estimates.) For all u\, u m , and u r G Bs 2 an d 1 < i,j ' < N 

the following property holds. Suppose that ui is connected to u m by an i-wave fan and that 
u m is connected to u r by a j-wave fan. Then, the wave strengths e k (ui, u r ) of the outgoing 
Riemann solution connecting Ui to u r satisfy (1 < k < N) 

e k (ui, u r ) = e k (u b u m ) + e k (u m , u r ) + 0(1) Q(u h u m , u r ), 

where Iguchi-LeFloch's interaction potential is defined by 

Q{ui 

) := ^Q^uuU^Ur) \o-i{ui 

i Urn) ®~ j{Umi ^r)| 

and 

{ 0, i < J, 



®ij{ui,u m ,u r ) :-- 



1, i > j, 

1, l = j, 0-j{u h U m ) &j(u m , U r ) < 0, 

Jj(u h u m ;u m ,u r ), i = j, a j (ui,Um)o- j (u m ,u r ) > 0. 



In the rest of this section, we derive a key estimate on wave speeds. Let us introduce 
the following notion. 

Definition 2.2. The inner speed variation of a j-wave connecting a left-hand state w„ 
to a right-hand state u + = ipj{m + \uJ), is defined as 

r dj{u_,u + ) = $j(m + ; uJ) 

t— max, . — min, , 

:= \j — A^- 

:= Xj(m + , m + ; u_) — Aj(m_, m + ; uJ). 

For instance, i?j-(it_, u + ) = if and only if w_ and u + are connected by a single admis- 
sible discontinuity, while ■#.,•(«_,«+) = \j(u + ) — Aj(w_) if (but not only if) w_ and u + are 
connected by a single rarefaction wave. 
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Theorem 2.2. (Properties of the inner speed variation.) 

1. For any u_, u'_, m it holds 

$j([j,(u'_) + m; u'_) = i?j(/i(tt_) + m; uJ) + 0(m) \u'_ — 

2. For any ui, u m , u r with u m = t/}j(^ m ] u{) and u r = ipj(fJ. r ] u m ), setting u' r := i])j{li r \ u{), 

$j{u h u' r ) = ^(/x P ;uj) = 'd j {n m ]Ui) + 0(\n r - /i m |). 

5. For any Ui,u m ,u r with u m = ipj(fi m ;ui) and u r = ipj(n r ]u m ) and monotonically 
ordered along the wave curve, 

l dj{Hr\Ui) <max u m ), $j(u m , u r )) 

( T 111 ' 11 / \ -r-min , , \ 

+ \\j {u m ,u r )~X j (ui,u m )J (2.2) 
+ 0(l)Q(u h u m ,u r ). 

In the application of Theorem 12.21 in the context of the front tracking scheme, every 
front («_, u + ) will propagate with the speed Aj" n (u_, u + ), so that the term (^\J m (u m , u r ) — 

\J m (ui,u m )j vanishes if the fronts (ui,u m ) and (u m ,u r ) interact. 

Finally, after establishing a density result of the class of non-degenerate flux (Theo- 
rem [231 below) we can conclude that: 

Corollary 2.1. All of the results stated in the present section remain valid for general, 
strictly hyperbolic flux that need not be nondegenerate, with the modification that a Riemann 
solution may contain (not finitely many but) countably many waves. 

2.3 Fundamental properties at critical points 

From now it will convenient to work with the following variant of the maps iTj, obtained 
by replacing the vectors r,j by the normalized vectors fj, i.e. 

9f(u) :=?,(«)■ VAjCu), 

?f +1 V) ;= f^) • VTrf (U), fc = l,2,..., 

Clearly, for every u and every j, the condition 

7rj 1) ( M ) = ... = 7rf( M ) = 0, 7rf +1) ( M )^0 

is equivalent to 

af («) = ■■■ = *f(«) = o, vf +1 Wo, 

so that a statement involving the coefficients ttj can be immediately restated with the 
coefficients 7Tj. 
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It is a well known fact that rarefaction curves and Hugoniot locus have a second- 
order tangency property at their base point. Here, we prove a new, higher-order tangency 
property which is satisfied at critical points. Recall that by differentiating the Rankine- 
Hugoniot it follows that 

(m — m_) d m Xj(m; w_) = Kj(m;u-) [Xj(vj(m; uJ)) — Aj(m;tt_)), (2.3) 

where the function kj = Kj(rn; u_) > is smooth, bounded, and bounded away from zero, 
and satisfies Kj(m_; u_) = 1. 

Lemma 2.1. (Tangency property at critical points (I).) Suppose that, at some point G 
Bg 1 and for some p > 1, 

= 0, k= l,...,p- 1, 

Then, the shock curve Vj and the rarefaction curve Wj issuing from u± are tangent at u± 
up to the order p + 1, 

dfchj(m*iu*) =d$w j (m i[ ;u i[ ), k = 1, . . . ,p + 1, (2.5) 

while the shock speed satisfies 

Proof. Proceeding by induction on we first note that, for all m, 

= Vnf\wj(m; w*)) • dmWjim; u*), (2.7) 

Now, (12.51) and (12.61) are proven by induction: we consider the induction hypothesis at 
the rank q (for < q < p): 

dm^Xjlm*, u*) = holds for k — 1, . . . , q — 1 and dm\v — w)(m*, u*) = holds 
for k — 0, . . . , q. 

This is clearly satisfied at rank q = 0. Suppose this is true at rank q < p — 1; let us 
establish it at rank g + 1. 

For the claim on dm Xj we proceed as follows. In view of (I2.3P and after differentiation 
with respect to m we obtain 

d^\ 3 (m- Ui ) + d$ ( ^'.^ SLA,-(m; = c^A^m; O). (2.8) 
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Considering (12.81) at the point m = m+, using (12. 7p and the induction hypothesis, it follows 
easily that dmXj(m+, u*) = for k = q with q < p. 

We now consider the second claim about dm +1 \v(m* , u*) — w(m* ,u*)) = 0. By differ- 
entiating the Rankine-Hugoniot relation, we find 

d m Xj(m) (vj(m) - uJ) = (A(vj(m)) - Xj(m)) d m Vj(m), 

which is to be compared with the equation characterizing the integral curve Wj(m) = 
Wj(m] it ), 

= (A(wj(m)) — Xj(wj(m))) d m Wj(m). 

The map z{m) := Vj(m; it_) — Wj(m\ u ) satisfies z(m ) = 0. From the above identities we 
deduce that 

(A( Wj ) - Xjiwj)) d m z ^ ^ 

= d m Xj (vj - u_) + (A(wj) - A(vj) - Xj(vjj) + Xj) d m Vj. 

Using (12.41) . (12. 7p . and the induction hypothesis on Xj (now proven up to the q-th order 
derivative) and by differentiating in m, the equation (I2.9p yields 

(A(u ) - Xj(u )) d^z{m*,u+) = 0. 
(In the case q = 0, no assumption on d$ Xj(m ir ,u ic ) is needed, but only Xj(m+;u+) = 



Xj(u+).) We deduce that, for some £ E, we have dm z(m*, u+) = 

Now, the parameter along the curves is such that 

fi(wj(m)) = m = fi(vj(m)). (2-10) 

Differentiating this relation and using, again, the induction hypothesis we obtain 

V/iK) d^ +1 S-K) = V^K) 0fi +1) Uj(m*), 

which yields dm +1 ^ z(m+, u+) = 0, and completes the induction. 

It remains to prove (12.61) at the rank p and (12.51) at the rank p + 1. Using Leibniz 
identity and retaining only the non-zero terms, we deduce from (12.81 ): 

d^XM; u0 + p d - X f m ^ = OtfXfaim*; «,)) = 7rj p) K), 
Kj[m, u+) 

which establishes the first claim. For the second claim, differentiating (12.91) p times, using 
Leibniz identity, and keeping the relevant terms only, we obtain 

(A(wj) - Xjiwj)) d^zim^u*) 

pd m Xj(m*,u±) - d m X j (w j (m^,u ir )) + ^Aj(m*, ii*)Jfj-(u*) = 0, 

in view of (12. 6p . Using again (I2.10p . this establishes (12.51) . □ 
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(2.121 



It follows from (I2.3P that, at a critical point m where 

a m A j (m ;n„) = 0, 
the shock speed must coincide with the characteristic speed, 

\j(m ;u-) = \j(vj(m ;uJ)), 
and it can also be checked from the Rankine-Hugoniot relation that 

d m Vj(m ;u_) = r,(^(m ; u_)). (2.11) 
Generalizing this observation we now prove: 

Lemma 2.2. (Tangency property at critical points (II).) If, for some «_,mo, and p, 

4 fc) Aj(m ;u-) = 0, k = l,...,p 

then at the critical point Uq := Vj(mo] uJ) one has 

7rf } (n ) = 0, fc = l,...,p-l, 
$rf («o)M 

(m - m_) ^ +1) Aj(m ; «_) = Kj(m ; u_) 7^- («o), (2.13) 

and £/ie Hugoniot curve issuing from u_ is tangent up to order p to the integral curve 
issuing from u , that is 

4^i(m ; u_) = uo), fc = 0, . . . ,p. (2.14) 

Proof. We first rely on ( 12.81) . By the assumption ( 12.12ft the left-hand side vanishes at 
m = m for all k < p. Hence, we deduce that, for all k < p, one has Xj(vj)(m ; u_) =0. 
With (jZIID we then deduce 

d m \j{vj{m Q ; «_)) = VAj(Uj-(mo; tt_)) <9 m ^(m ; u_) = VAj • ^(^(mo; u_)), 

Hence, (i>j-(mo;tt-)) = 0. Hence, using ( 12. 9ft and ( 12.101) . one sees, as in Lemma [27TI that 
(l2TT2j) is valid for p = 2. 

More generally, if (I2.14p is already established at the rank k then using ( 12.71) it follows 
that 

d m x j ( v j (™o! «-)) = ^S? Aj (wj (m ; M ) ) 

Hence, we conclude that 7rjj (ito) = 0. Using the same procedure as above, one deduces 
(J2H3D at the rank k + 1. 

For k = p the left-hand side of (12.81) at m is a multiple of the left-hand side of (12.131) . 

m .°- m - e i} AiK;u-), 

while the right-hand side equals 7T^(m ). Thus, (12.131) holds. □ 
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Some further observations are in order. 

Lemma 2.3. (Propagating discontinuities with coinciding speeds.) Suppose that (uq,ui) 
and (ux,u 2 ) are shock waves satisfying the entropy criterion and propagating with the same 
speed A, 

Ul = Vj(mi;u ), U2 = Vj(m 2 ;ui), 

\j(u , ui) = Xj(ui, u 2 ) = A, 

where mi and m2 satisfy h(uq) < ni\ < m 2 . Then, the discontinuity (uo,u 2 ) is a shock 
satisfying the entropy criterion and propagating with the speed A, 

u 2 = v j (m 2 ;u Q ), X j (u ,u 2 ) = A. 

Lemma 2.4. (Equivalent formulation of the entropy criterion.) A discontinuity («_,«+) 
with u + G Ti.j(u-) satisfies the entropy criterion if and only if 

Xj(u-, u + ) < Xj(m; U-) for all m between m_ and m + , 

where m_ = /i(uJ) and m + = (J,(u+). 

We only give the proof of Lemma 12.41 the proof of Lemma 12.31 being similar. 

Proof. Set m_ := ji{u_) and m + := yu(u+) and, for definiteness, assume that m_ < m + . 
We must show that 

A := Xj(u-,u+) > Xj(m;u + ), m_ < m < m + . 

By contradiction let us assume that there exists a "first" point m G [m_,m + ) at which 
the above condition fails, that is 

A > Xj(m; u + ), m_ < m < mo, 

A = Xj(m;u + ), m = mo, (2-15) 
A < Xj(m; u + ), < m - m <C 1, 

On the other hand, since the shock connecting U- to u + being entropy admissible, 

A < Xj(m;uJ), m_ < m < m + . (2-16) 

We treat the case where m_ < mo < m + , the case where the point mo coincides with the 
left endpoint m_ of the interval being analogous. 

Let us expand Xj(m; u + ) in a neighborhood of m = m , 

Xjim; u + ) = A + ^ ~ p ™^ ^ +1) A,(m ; u + ) +0((m- m ) p+2 ) , 

where dm +1 ^ Xj(m ; u + ) ^ 0. By Lemma l2~2l we know that up to a positive multiplicative 
constant 

dt +l) \ J {mo;u + )=9 i f\uo), 
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while all 7?j (uq) = for k = 1, . . . ,p — 1. On the other hand, since the shocks (u_, u + ) and 
(u , u + ) propagate at the same speed, A, the shock speed of uq) is also A. (This is an 
elementary fact, already stated in Lemma [2.31 ) Therefore Aj(m ;u_) = A and, similarly 
as above, we can write 

Xjim; u_) = A + {m ~^* 1 e^Vmo; «-) + 0((m - m )«+ 2 ), 

where 

^ +1 >A,(m ; M _)=^W 

while all 7?^ (wo) = for A; = 1, . . . , q — 1. Obviously, q — p. However, in view of conditions 

(12.151) on Xj(m]u + ) we see that p must be even (with ttj (uq) > 0), while the condition 

(12. 161) on Xj(m; u_) implies that p must be odd (with 7r) | (mo) > 0). This is a contradiction. 

□ 



2.4 Construction of the wave curves 

In the rest of this section we sketch the proof of Theorem 12. 11 Most of the arguments in [23] 
carry over under our weaker assumption on /, and we will not repeat them. We will only 
give the new ingredients of the proof. For defmiteness, we assume that (VXj ■ fj)(u ) > 0, 
and we construct the part m > of the wave curve. Locally near u we can use the 

integral curve Oj(uo). 

Suppose that along the integral curve there exists a "first point" with coordinate ^(uq) 
where VXj ■ Tj vanishes and changes sign, that is by setting u 1 := Wj (fj} (u ) ; u ) 

(S7Xj ■ r^iwjim; u )) > 0, fi(u ) < m < ^(uq), 
(VA r f J -)(« 1 ) = 0, 

(2 17) 

(VAj •rj){wj{m;uo)) < 0, < m - fi (u ) « 1, 
7^V) = 0(* = 1,...,P), ^ +1) (n 1 )<0. 

Note that, clearly, p must be an odd integer. 

It should be noted here that there can not be accumulation points in the critical set 
along an integral curve. Indeed, if m —>■ m°° (I — > oo) were a sequence such that 

(VXj ■ f;-)(w,-(™ l ;«o)) = 0. 

Then, by induction on k > 1 and using the intermediate value theorem we could find 
sequences m k ' 1 such that 

Letting / — > oo we would conclude that, at the point u°° := Wj(m°°; Uo), 

rrf\u°°)=0, fc = l,2,..., 



11 



which is impossible since / is nondegenerate. 

Note that the wave curve coincides with the integral curve up to the value /x 1 (w ): 

ipj{m; u ) := Wj(m; u ), fJ,{u ) < m < ^(uq). 

We claim that fi 1 is a continuous function of its argument. Namely, consider the function 

Fj(m; u) := 7rf\wj(m; u)), 

which is smooth with respect to m and u. Since by assumption the state Uq is such that 
FjQi?- (uq);uo) = 0, 

(d^F 3 )(^(u );u ) = ^ +1) { Wj {m-u)) < 0. 

Therefore, one can locally solve the equation Fj(fi,u) = in a neighborhood of uq, in a 
continuous manner. 

To extend the wave curve we need a pattern made of a rarefaction wave followed by a 
shock wave (actually a left-contact wave). For all meaningful values (n, m;u) we set 

(Aj(m; Wj(n; u)) — \j(wj(n; u))) , m ^ n, 



Gj(m,n;u) := <{ f~ n 



2 

where u 1 := Wj(/j}(u);u). We will now apply the implicit function theorem to the equation 

Gj(m,n;u) = 0, (2.18) 

and show that it defines a function n = z/ 1 (m; u) in the neighborhood of the point (m, u) = 
(/x 1 («o), uq). As presented above, our construction is associated with the point u 1 at which 
r j ' VA^u 1 ) = 0. The special case of a PGNL flux is straightforward since it can be checked 
that rj ■ V- (?".,• ■ VAjXu 1 ) ^ implies that 9 n G J (/i 1 (u ), /i 1 (u ), u ) ^ 0. Handling a general 
nondegenerate flux is much more delicate and this is discussed now. We regard u as a 
parameter and, for convenience, in the presentation we simply fix the point uq and skip 
the dependence in u. We are interested in the corresponding singular point u 1 along the 
rarefaction curve from w , where n^ p+l \u l ) ^ and we must analyze (12.181) near 

m 1 ■= ^(uq). 

Consider first the case of a scalar conservation law. The function Gj is independent of 
the variable u, and the equation under consideration reads 

" /M - f\n) = 0. 
n — m 

Using that / has a nondegenerate point at m 1 , precisely, f^ p+2 \m}) ^ so that 
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one can check geometrically on the graph of the function / that the above condition 
determines a unique solution n = i/ 1 (m) which, furthermore, can be easily expanded near 



m 1 : 



z/ 1 (m) ~ c p [m — m 1 ] 



where the constant c p < is the unique negative root of the equation 

1 - c p+2 

-—L- = {p + 2)d*\ 

1 Cp 

We now extend the above result to systems. 
Lemma 2.5. For a, (3 > sufficiently small, introduce the small cone 

C := ||(n — m 1 ) — c p (m — m 1 )] < a \m — m 1 )}. 
Near the point (m, n) = (m ,m ) within the cone C the function Gj admits the expansion 
Gj(m,n) = Q^m.n) + 0(1) ((m - m 1 )^ 1 + (n - m 1 )^ 1 ), (2.19) 

where 



(2.20) 



Qj(m, n) = (jji — m l ) — c p {m — m 1 )) Rj(n, m), 
Rj(m, n) = ^ +1 \u x ) (m - m 1 )^ 1 Sj(n, m), 
and the function Sj is smooth and close to 1, 

\Sj{m,n) - 1| < p. (2.21) 

Once this lemma is established, the implicit function theorem straightforwardly applies 
to the mapping Gj defined by 

— (jr ■ ( m n ) 

Gj{m, n) := 3 [ in the cone C, 

{m — m 1 -) p 1 bj{m,n) 

and extended arbitrarily outside C as a smooth map. By the above claim, we have 

Gj(m,n) = 7T (p+1) (w 1 ) ((n - m 1 ) -c p (m- m 1 )) + 0(l)(\m - m 1 ] 2 + \n- m 1 ] 2 ), 

so that there exists a unique function u 1 = i/ 1 (m) defined in a neighborhood of m = m 1 
and such that 

Gj{m, v 1 {m)) = 0, 

which, moreover, satisfies 

^ 1 (m) = m 1 + c p (m — m 1 ) + o(m — m 1 ). 

In turn, the solution remains within the cone C (in a sufficiently small neighborhood of 
m 1 , at least), and therefore i/ 1 (m) is a solution of the original equation f)2.18p . 
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It remains to prove Lemma 12.51 above. Relying on the tangency property between the 
Hugoniot and integral curves (Lemma 12.11) . one can replace the Hugoniot curve by the 
rarefaction curve in the definition of Gj while making an error of order (m — m 1 ) p+1 + (n — 
m 1 ) p+1 . Let us introduce the reduced flux 

fj(u) :=tf(u) 

and set 

%{m) := fj{wj{m;uQ)). 
We can rewrite the equation Gj(m,n) = in the following form 

fj{m)-fM) = ju n) + ,y n m _ m iy+i + r n _ m y+i). 
m — n J 

The behavior of the flux in the neighborhood of u 1 leads to 

and the expression in (12.171) can be expanded exactly as in the scalar case, leading to the 
result stated in the claim. 

The mixed curve is now defined locally, and we discuss its extension. Consider the 
shock speed of an arbitrary discontinuity connecting Wj(n; u ) to Vj(m; Wj(n; u )), that is, 

Aj(m, n) := Xj(m;wj(n;uo)). 

By construction, at n — z/ x (m) the shock speed coincides with the characteristic speed: 

Aj(m, v l {m)) = Xj(wj(iy 1 (m);u )). 

One can show that 

d m u\m) = — — — -. 1 rr (d m Aj)(m, v\m)). 

(VAj •r i )(w i (z/ 1 (m);M )) 

Since v l \[i l ) < 0, it follows that v v (m) < for m sufficiently close to /i 1 and therefore 

(d m Aj)(m, //(m)) < 0, < m - fj}(uo) « 1. 

This means that the speed of the left-contact decreases as m increases. The function 
i/ 1 (m;w ) is well-defined until d m Aj eventually vanishes (see in particular d m Gi) at some 
point with coordinate denoted by fi 2 = /i 2 (w ), 

(^AjX/xVV 2 ;^)) =0. 
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The wave curve is determined from the map m Vj(m; w : ,(i/ 1 (m; u ); u )), by 

\wj(m;u ), fji(uo) < m < ^(uq), 

As in [23] one can check that for each m G (jt/ 1 , fi 2 ) the state w J (z/ 1 (m); «o) is connected to 
the right-hand state Vj(m; iw^i/^m); uq)) by a shock wave satisfying the entropy criterion. 
All of the remaining arguments of the construction of the wave curve in [23] can be extended 
in a similar way as we have explained above. 

We observe that the construction requires finitely many waves only. Otherwise, if the 
wave curve was made of infinitely many shock and rarefaction pieces, then there would be 
an accumulation point of critical points, m l — > m°° (I — > oo) such that 

7rj 1) (^-(m I ;«o)) = 0. 

Then, by induction on k > 1, using the intermediate value theorem and the higher-order 
tangency property between shock and rarefaction curves as stated in Lemma 12.21 we would 
obtain 

for all k, which is impossible since / is nondegenerate. 

After the wave curves are constructed, the Riemann problem can be solved as in [23] . 
This complete the discussion of Theorem 12. II Propositions 12. II and 12.21 are then immediate 
from the results in [23]- From our construction, the following property of the Riemann 
solution follows immediately. 

Lemma 2.6. (Splitting property for the Riemann problem.) Consider a wave fan u + = 
ipj(m+;uJ), with rri- < m + for definiteness, associated with the intermediate states Uq, 
. . . , U2 q +i as defined in Theorem \2.1[ Then the following two properties hold. 

1. For k = 0,...,q and m G [/x(ii 2 fc), M w 2fc+i)]> ^ e solution of the Riemann prob- 
lem (u_,tpj(m;u_)) (resp. (tpj(m;u_),u + )) coincides with the restriction of the solution 
of the Riemann problem (u-,u+) to the interval [Aj(m_, m + ; \j(i/jj(m; U-))] (resp. 
[\j(ipj(m; it_)), Xj(m + , m + ; uJ)\). 

2. For k = 0, . . . , q and m G [fi(u2k+i), ^(^2^+2)]; the solutions of the Riemann prob- 
lems (u-,i/)j(m]U-)) (respectively (^(m; w_), u + )) and (m_,m + ) coincide in the interval 
[Aj(m_, m + ; u_), \j(u 2k+1 )] (respectively [Xj(u^ k+2 ),Xj{m + , m+; U-)\). 

We end this section with another property of Riemann solutions: roughly speaking, the 
speeds within a wave fan can be determined by a convex or concave hull argument, which 
is similar to what is classically done for a scalar conservation law. This property will not 
be directly used in the rest of this paper. 

Lemma 2.7. (Construction of a wave packet using the convex hull.) Let u + = ipj(m + ; U-), 
with m_ := yu(u_) < m + , and introduce the intermediate states U\,. . ■ ,u 2q +i as in Theo- 
rem \2.1\ By setting 

fj(m) := I ■ f(ijjj(m] u_)), m_ < m < m + , 
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the convex hull gj := conv[ m _ >m+ ]fj is precisely 

fj( m ) = fj( m 2k) + / \j(ipj(m'; u_)) dm', 

Jrn 2k 

gj (m) = ^ m 2k <m< m 2k +i, (2.22) 

fj(m 2 k+i) + Aj(« 2 fc+i,«2Jfc+2) (m - m 2 jfc+i), 

m2t+i < m < m 2 fc+2- 

Moreover, fj cannot be affine in a suhinterval o/[m_, m + ] where it coincides with conv\ m _ )m+ ]fj 
and, in consequence, the set of the intermediate values m k is nothing but the boundary of 
the set 

{m e [m-,m+]/ gj(m) = fj(m)}. 

When m + < m_, the same statements hold by replacing the convex hull by the concave 
one. 

Note that, even though the wave curves ipj are only Lipschitz continuous, the function 
gj above is of class C 1 , as this is clear from its definition. 

Proof. Let us here denote by gj the right-hand side of f)2.22p (defined successively on 
each interval [m k ,muk+i\- We will show that gj = convr m _ )m+ ] fj. Observe that gj is a 
convex function, as a consequence of the fact that the wave speed function \j(m, m + ; uJ) 
introduced earlier is non-increasing in m. 

First of all, within a rarefaction interval [ffi2)t) w*2i+i] t ne identity 



fj(ipj(m,u-)) = fj(m 2k ) + \j(ipj(m';u-))dm' (2.23) 



follows directly by differentiating this relation and using the definition of the rarefaction 
together with the normalization fi(vj(m)) = m. 

Second, considering a shock wave (u 2k , u 2k+ i) and using the Rankine-Hugoniot relation, 
we obtain (for the end point values) 

fj{u 2k+2 ) - fj(u 2k+1 ) = Xj(u 2k+1 ,u 2k+2 )T- (u 2k+2 - u 2k+ i) 

= Xj(u 2k+1 ,u 2k+2 ) [m 2k+2 - m 2k+l ). 

Hence, combining with the identity within rarefaction, we conclude that fj and gj coincide 
within rarefaction intervals, 

fj(m) = gj(m), m e [m 2k ,m 2k+1 ]. (2.24) 

In consequence, to establish (12. 22ft it is sufficient to prove that gj < fj. Moreover, 
since these functions coincide on the intervals [m 2k , m 2k+ i], we need only show gj < fj 
on each interval [m 2k+ i,m 2k+2 \. So, we now establish that, for any u = i/jj(rh;u_) with 
m 2k+ i <rh< m 2k+2 , 

fj(m) > fj{u 2k+ i) + \j(u 2k+1 ,u 2k+2 )(rh - m 2k+l ) 

/~x {2. 2d) 

= 9j{m). 
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To this aim, we consider the Riemann problem («_,«) and the restriction of fj to the 
interval [m_ , fh] , 

fj - = fj\[m-,fh] 

and we define fjj by (12.221) with the intermediate states u[, . . . , u' 2 , +1 corresponding to this 
problem. It follows from previous observations that, at the end point m, 

fj(rh) =gj(m). 

In view of the splitting property in Lemma |2~B"1 it is clear that gj(m 2 fc + i) = fj(m 2 k+i), an d 

9j(rn) > fl5(m 2 jfe+i), m 2k+l <m<fh. (2.26) 

Now, if u 2 k+i 7^ u~, we have 

9j(m 2k+1 ) = fj(u 2 k+i) 

= \j(u 2 k+l) = \j{U2k+l, u 2k+2), 

since the two Riemann solutions coincide on a whole interval. If u 2 k+i — we have 

gj{™>-2k+l) = ^j(u 2 k+l,u) >\j{u 2k+ i,U 2k+2 ), 

as follows from the entropy criterion. This establishes (12.251) . 

Consider now the last statement in the lemma and suppose that fj is affine and co- 
incides with conv/j in some subinterval [mi, m 2 ] (mi < m 2 ). Clearly, (m 1 ,m 2 ) does not 
intersect a rarefaction interval (m 2 k,m 2 k+i) since this would contradict the nondegener- 
acy assumption. Introducing u as above, we see that the inequality (12.261) is strict unless 
(«2fc+ij u) is a single shock of the same speed as (u 2 k+i, u 2 k+ 2 ). But, if Vj(-; u 2 k+i) had this 
property on [mi,m 2 ], it would follow that Vj(-; u 2 k+i) = Wj{- \ u 2k +i) in that interval which, 
again, would contradict the nondegeneracy of the flux. □ 

2.5 Proof of the inner speed variation estimates 

Proof of Theorem 12.21 1. We only treat the case m > since the case m < is 
similar. From the estimates on the strengths we have: 

sup |'0 J -(//(it_) + n; iz_) — tl)j{^{u'_) + n; u'_)\ < |m| \u'_ — (2.27) 

0<n<m 

Observe also that the variation of the shock speeds is essentially equiavelent to the variation 
of the characteristics speeds, as follows 

Aj'OOO + n; u') - \j(n(u) + n- u) 
= Aj(m') - Xj(u) + 0{n) \u' - u\. 



20 



This is a consequence of the fact that all functions under consideration are smooth, and 
that the equality holds (without remainder) when either n = or u' = u. 

Denote u' + := ipj(fi(u'_) + m;u'_) and u + := ?/^(/i(u_) + m;w_). To prove the first 
statement in Theorem 12.21 we will establish 

A™ aX (/x(w'_) + m; u'J) — A™ ax (/i(n_) + m; u„) 
= \j{u' + ) — \j(u + ) + 0{m) \u' + — u+\ 



(2.29) 
(2.30) 



{fi{u_) + m;u_) - Aj {fi{U-) + m; tt_J 
= Aj(ii'_) - Aj(w_) + O(m) - 
This implies the desired estimate, since by (12.271) . 

Aj(«+) — Aj(w + ) = \j(u'_) — Xj(u-) + C(m) — U-\ 

and — u+| = \u'_ — u~ \ + 0(m) — u_|. 

We only consider (12.291) . since (12.30jl can be checked similarly. To prove ( 12.291) . it is 
sufficient to establish the inequality 

\J ax (n(u'_) + m; «'_) - A7 ax (/i(«_) + m ; «-) ^ ^ 

< Aj(w' + ) — Aj(m + ) + 0{m) \u' + — u + \. 

This is so since u_ and u'_ play completely symmetric roles. 

Each of the Riemann problems («_,«+) and (u'_,u' + ) are solved by a succession of 
rarefaction waves and shock waves. Call m (resp. m') the real in [0, m] defined as follows: 

• If the problem (u_, u + ) (resp. (u'_, u' + )) ends with a non-trivial rarefaction wave, we 
set m = m (resp. m' = m), 

• If the problem (u_,u + ) (resp. (u'_,u' + )) ends with a non-trivial shock wave (u, u + ) 
(resp. (u, u' + )), we set m = ji{u) — ji{u~) (resp. m' = /i(£t) — /i(n'~)). 

We distinguish between two main cases: 

• Case 1: if m' > m. By the entropy criterion, the shock speed enjoys a monotonicity 
property and we have 

Aj(/x(m_) + m; u + ) < Aj(//(«_) + m; u+) = A™ ax (n_, 

and the conclusion follows easily from (12.281) . 

• Case 2: if m' < m. In that case, there are two subcases: 

— 2a. If the point u of parameter value /x(tt_) + m' in the problem («_,«+) 
corresponds to a rarefaction wave. In that case, we clearly have 

Xj(fi(u_) + m, yu(n_) + m; u_) > Aj(w). 

Hence, using Lemma [2.61 we can prove ( 12.31}) as in Case 1, by considering the 
Riemann problems (it_, if)j([j,(u-) + m/;u-)) and {u'_,ij)j{^{u'_) + rn';u'_)). (In 
fact, in this case we only compare characteristic speeds.) 



21 



— If the point u of parameter value //(«_) +m' in the problem (w_, u + ) corresponds 
to a shock wave, say (ui,w 2 ), then 

Aj(/x(tt_) + m, + m; tt_) > A J (ti 2 ). 

Again, using Lemma 12.6} we can apply Case 1 to the problems (tt_,tt2) and 
(?/_, tpj(fi(u'_) +fi(u 2 ) — «'_)). Note here that the latter problem coincides 

with (u'_, u' + ) up to the parameter fi{u'_) + m' and has speeds greater than the 
ones in (u'_, u' + ). 

2. Fix a left-hand state Ui, and let u m , u' r be two states on the j-wave curve from ui, with 
u' r := ipj{fi r ] ui). Using the standard Lipschitz continuity property of the wave speed along 
a given wave curve based at the point ui we obtain 

A™ ax (wi, u' r ) = Xj ax (u h U m ) + C(/i m - fl r ), 

-r-min , . . -r-min , . . . 

Aj [Ul, U r ) = Xj [Ui, U m ) + U{fl m -fl r ), 

which yields the desired statement on the inner speed variation $j(ui, u r ) = A™ ax (n;, u' r ) — 
Af n K<). 

3. Introduce u r '■= ipj{fi r '-, ui)- As we could construct "left-hand" wave curves as we did for 
(right-hand) wave curves, one can introduce u m such that 

It follows from the regularity of the left-hand wave curves and the wave interaction esti- 
mates that 

\u m - u m \ = 0{\u r - u r \) = 0(l)Q(u h u m ,u r ). 
Now, using the first statement of the theorem, we deduce 

Xj(-, fir] U m ) = Xj(-, fl r ] U m ) + 0(1) Q(u h U m , U r ). 

We introduce two values of the parameter in the Riemann problem (ui,u r ): 

m a := maxjm G [fii, fi m ], / the point in the problem (ui,u r ) 

with parameter m corresponds to a rarefaction wave}, 

Tiib := min \m G [/i m , /i r ], / the point in the problem (ui, u r ) 

with parameter m corresponds to a rarefaction wave}, 

with the convention that m a = \xi (resp. m& = fi r ) when there is no rarefaction wave in the 
corresponding parameter range. In the terminology of Lemma 12.71 the above conditions 
can be stated as fjijn) = conv[ W)Atr ]/j-(m). 

We now distinguish several cases that depend upon the values of m a and ra&. 



22 



• Suppose that \i\ < m a < m h < \i r . In that case, using Lemma [231 we see that the 
Riemann problem (ui,u r ) can be obtained by gluing three Riemann problems to- 
gether: (ui,if)j(m a ;ui))i (ipj( m a',ui),ipj(mb;ui)) (a single contact discontinuity), and 
(ipj(mb,ui),u r ). We then deduce that 

■&j(ui, U r ) 

= Xj(fl r , Ur] U m ) - Hm\ Ui) 

— Xj{Hrj Hr\ Urn) — Aj(/i m , \l r \ U m ) + Xj(fJ, m , fl r ] U m ) — Xj(^l, Hm'i Ul) 

my f^rj ^m ) - fi m ; ui) + 0(l)Q(ui 

and the desired inequality ( 12. 2ft follows. 

• Suppose Hi = m a < m b < /i r (resp. Hi < m a < m b = Mr)- I n that case, by 
the same arguments, we have $j(ui,u r ) = , dj(^j(mb;ui),u r ) (resp. $j(ui,u r ) = 
r dj{ui ) %ljj{m a ]Ui))) and, again, ( 12.21) follows. 

• Suppose hi = m a < mb = fir- This case is obvious since the left-hand side of ( 12. 2 j) 
vanishes. 

This completes the proof of Theorem 12.21 □ 



2.6 Proof of the density property 

We will now establish that the nondegeneracy condition is generic. This section will cover 
both conservative and nonconservative systems, characterized by a flux / or a matrix A, 
respectively. 

Let fl be an open set in R N . We are interested in matrix-valued mappings A = A(u) 
(u G fi) of class C°° satisfying the strict hyperbolicity property, that is, for any u G Q, the 
matrix A(u) admits N distinct and real eigenvalues, 

Xi(A,u) < ... < X N (A,u), (2.32) 

and basis of left- and right-eigenvectors lj(A, u), rj(A, u), normalized so that 

l r i(A u)\ = 1, lj(A,u) -rj(A,u) = 1, 
k(A J u)-r j (A J u) = (i^j). 

Let H, be the subset of C°°{Vt) consisting of all mappings A satisfying the strict hyperbolic- 
ity condition (I2.32p and such that the associated maps Xj, lj,rj are smooth. Recall that, in 
a sufficiently small neighborhood of any constant matrix and under the assumption (12.321) , 
the maps Xj,lj,rj indeed depend smoothly upon u. In this subsection however, we need 
not assume that the set Q is small. 
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To every map A G 7i we associate the following functions, by induction, 

(2.34) 



irf\A, u) := rj (A, u) ■ V u \j{A, u), 



vrf +1 \A,u) := rj (A,u) ■ V u nf\A,u), k = 1,2, . . . 
We claim that: 

Theorem 2.3. (Density of the set of nondegenerate systems.) Given any matrix-valued 
map A G Ti there is a sequence A 1 G 7i such that 

A 1 — > A in the (strong) C°° Whitney topology, (2.35) 

and A 1 is nondegenerate in the sense that, for I — 1, 2, . . ., u G tt, and j = 1, . . . , N, 

(nf(A l , u), Trf (A 1 , «),..., vrf +1) {A\ «)) ^ (0, 0, . . . , 0). (2.36) 

Moreover, in the special case that A is conservative, i.e. A = Df, the sequence can be 
chosen to be conservative, A 1 = Df 1 , and f l converges to f in the Whitney topology. 

Hence, all the results established earlier in this section for nondegenerate flux extend 
to arbitrary, strictly hyperbolic flux, and this leads us to Corollary 12. 11 In the following, a 
map A G Ti (conservative or not) is fixed, and we consider the set 

T := {(u,A(u)), u G tt}. 

Clearly, (12.321) is an open condition and, therefore, for each u G tt one can find S > such 
that the ball Bs(u) centered at u is included in tt and all matrices B G B$(A(u)) have N 
real and distinct eigenvalues: 

A 1 (B)<...<A N (B), 

and, therefore, basis of left- and right-eigenvectors Li(B), . . . , L N (B) and Ri(B), . . . , R N (B), 
which we will normalize so that 

\Rj(B)\ = l, Li(B) ■ Ri(B) = 1, 

Li(B) ■ Rj(B) = (i^j). 1 • ) 

Clearly, the maps Aj,Lj,Rj depend smoothly upon B (at least locally), and represent 
extension of the maps Xj,lj,rj to the ball B$(A(u)). Reducing 5 if necessary, we can 
assume that for all Bi, B 2 G Bs(A{u)) 

Li{B x ) ■ Ri(B 2 ) > 1/2. 

Now, in view of the paracompactness property of J 7 , there exists a locally finite covering 
of T by open sets finer than Bs/2((u, A(u))). We define Q as the unions of all balls in 
such a covering. We define Rj(B) and Lj(B) on Q for B G Bs/2(A) as the value of 
Rj{B) determined in that ball. It follows easily from this construction that this does not 
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depend on the ball that contains B. Hence, we conclude that in Q the maps Q 3 A i— > 
Aj(A), Lj(A), Rj(A) are smooth. 

To fix the idea we discuss the conservative case, with A = Df for some flux /. Minor 
modifications are needed to cover the nonconservative case, which we omit. 

Let O C C°°(Q, WL N ) be a small neighborhood of the flux /, chosen so that (u, Dg(u)) £ 
Q for all g £ O and u £ Q. Introduce the subset 

J k (0) = {j k g (u) := (u,g(u),g^(u),...,g^(u)) / u e tt, g e O} 
C J k (Q,R N ), 

consisting of all k-th jets of maps from into M. N . Here, g&> denotes the j-order differential 
of the map g, and the mapping j k g is the k-th jet extension of g. 

It follows easily from the definition (12.341) that the functions u i— > Kj (/, u) can be 
expressed in terms of jets, that is 

*f (/, «) = /(«), • • • , (2 3g) 

= Vfc(jJ +1 («)), 

where the functions (p k '■ J h+1 (Q,R N ) — > K are smooth in J7a;+i(0). Note in passing that 
ftj k \f, u) does not depend upon the first two components of the jet, u and f(u). 

The proof of Theorem 12.31 follows from Thorn's transversality theorem, which we recall 
for readers' convenience (for a proof see, for instance, [22]), and a technical lemma that we 
establish below. 

Theorem 2.4 (Thorn). Let X and Y be two smooth manifolds, and Z be a submanifold of 
the k-th jet space J k (X, Y) for some k > 0. Then, there exists a Baire set of second category 
E C C°°(X, Y) for the (strong) C°° Whitney topology such that the k-th jet extension of 
any f £ E is transverse to Z, and consequently (jj) _1 (Z) is either the empty set or a 
submanifold of M having the same codimension as Z. 

Recall that a Baire set is a set whose complement is the union of at most countably 
many, nowhere dense sets. 

Lemma 2.8. For any k > 2, the set 

Z:=[qe J k (0) I tpfa) = 0, j = 1, . . . , k - l} 

c j k {Q,m N ) 

is a submanifold of codimension k — 1 . 

Proof of Theorem 12.31 Applying Lemma [2781 with k = N + 2, we obtain a submanifold Z 
of J N+2 (Q, R ). By Theorem 12.41 the set E of maps g : Q — ► W N such that g is transversal 
to Z is a Baire set of second category. For g £ E we introduce 

^ g :u^(nf\g,u),...^ k - 1 \g,u)). 
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For any g £ E, 1 (0) is either empty or is a submanifold of codimension N + 1. Clearly, 
there are no submanifold of codimension TV + 1 in Q C R , hence \E f ~ 1 (0) = 0. Thus, one 
can find a sequence / G E converging to / for the Whitney topology. □ 

Proof of Lemma 12.81 It is enough to prove that the differential dip k of the mapping 

ij k : J k (Q, R N )3q^ ( Vl (q), tp h -. x {q)) 

has constant maximal rank k — 1 for each q in J k (tt,R N ). Let us write (A , A u . . . , 
for a general vector of the tangent space at q, T q J k (Q, R N ), with A ,Ai G R N , A^ G 
i > 1, where dj = (JV + i - 1)\/[(N - 1)H\}. 

By a straightforward induction argument we can check that it is sufficient to prove, for 
every k, 

(0,...,0,1) G Range(DVfc) ■ (2.39) 
Denoting by the fc-th order differential of /, we first observe that 

^(f,u) = D u [A i (f^(u))}(R i (f( 1 \u))) 

= (D A Ai)(fW(u)) ■ (/ (2) )0) • [Ri(f {1) (u))], 

where p 2 > is regarded as an element of C(M. N ; £(R N ,R N )). 

It is straightforward to check, by induction, that the term in (f k -i containing the 
highest-order differential of / is 

(£> A A<)(/ {1> («)) • / (fc) («) ■ [^-'Riif^iu))}, 
where f^ k > is regarded as an element of C((R N )® k ~ 1 ; C(R N ,R N )). In other words, we have 

V7fc-l(«, /(«),..., /«)(«)) 

= (D A A i )(f^ 1 \u))^ k \u).[^ k - 1 R i (f^(u))} + ^(u, /(«), . . . , Z^- 1 )), 

where ^fc-i is a smooth function of 

Hence, to establish (12.391) . it is enough to prove that for all A G 7i there exists 

A k G £((R N f k - l ;£(R N ,R N ))^ £((R N f k ;R N ) 

such that 

D A Ai(A) ■ A k ■ (^ k - l Ri(A)) 0. (2.40) 
Note that A, Ri(A) and Ai(A) are related together via the identity 

ARi(A) = Ai(A)Ri(A). 

Differentiating this identity with respect to A (in the direction if), we obtain 

H Ri(A) + A- D A Ri(A) ■ H 
= (D A A t (A) ■ H) Ri(A) + A t (A) D A R t (A) ■ H. 
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Taking the product by the left eigenvector Li(A) and using (12.331) . this yields 

D A A l (A)H = L t (A)HR t (A). 

Hence, we arrive at (12.401) by considering A^ sending ® k Ri(A) to Ri(A) and other elements 
of a basis of (R N )® k to 0, and A = 0,. . . , A k . x = 0. □ 

Remark 2.1. It follows from the above proof, that generically in f (or in A in the non- 
conservative case), one has the stronger property that 

{uen/ (nf(u),...,nf\u))} 

is a (possibly empty) submanifold of codimension k. 

As a consequence, for such a flux and for almost every u£fl, the (Lipschitz continuous) 
wave curve starting from u does not meet any point where (tt^\u), 7rj (w)) = (0,0). 



2.7 The PGNL condition is not generic 

We have established that any flux can be approached by a sequence of flux satisfying 



0) 



for all u and j. We will now prove that, in general, a flux cannot be approached by a 
sequence of flux satisfying the stronger condition 



(u))^(0,...,0) 



for all u and j. Hence, the latter condition is not generic. In particular, for N > 2, there 
exist fluxes that cannot be smoothly approached by PGNL fluxes. Note that our coun- 
terexample below is "local" , in the sense that it persists even if we shrink the neighborhood 
of the base point under consideration. Moreover, our result is valid for both conservative 
or nonconservative hyperbolic systems. 

Consider flux / : M. N — > ~R N of the following triangular form: 



/ 



whose Jacobian matrix A(u±, . 



,u N ) 



/i(«i) 



\ 



( 



d Ul f2(ui,u 2 ) 



\ /jv(«i, • 
,u N ) is given by 


^2/2(^1,1*2) 



wjv) / 








\ d Ul f N (u!,...,u N ) d U2 f N (ui, . . . ,u N ) 



d UN fN{ui, . . . ,u N ) J 



27 



Under the assumption 

duifi < 9 U2 f 2 < . . . < d UN f N , 
the system of conservation laws associated with / is strictly hyperbolic, with 

/0\ 

Xj(u 1 ,...,u N ) = d Uj f j (u 1 ,...,u N ), r N (u 1 ,...,u N )= ■ 

V 1 / 

It will be convenient to choose the component to be the following polynomial ex- 
pression 

f N (ui, u N ) := Mi u 2 N + u 2 u N + . . . + u N -i Mjy + m^ +2 , 

while fi, . . . , /jv-i are chosen arbitrarily but so that the strict hyperbolicity property holds 
for all u in the ball B$. One can compute 



Ajv(ui, 



MAT 



= 2 Ul u N + 3u 2 u^ + . . . + Mijv-iujJ -1 + {N + 2)«jJ +1 , 
thus, since • V = d UN , 

TTx'{ux, ■ ■ -,u N ) := r N ■ VAjv(wi, . . . ,itjv) 
= 2ui + . . . + N(N - 1) Mjv _i <" 2 + (N + 2)(N + 1) u#, 



(2.41) 



(2.42) 



and so on for all nj , in particular 



7rjf 1} (tti, • • • ,u N ) = N\ u N ^ + 
7$° ujv) = (iV + 2)! n,v. 



(JV+2)! ,9 



AT) 



(2.43) 



Now, we claim that: 



Proposition 2.3. Consider the triangular flux f introduced above. For any 5 G (0, 1) and 

for all g sufficiently close to f on Bs (in the C N+1 -norm), there exists a state u* G B$ such 
that 

(v4V),---,^K)) = (0,...,0). (2.44) 

Our proof below remains valid if, instead of flux g, we consider mappings A = A(u) 
(u G Bs) which are sufficiently close to df. 



Proof . Clearly, any perturbation g : B§ 



of / is still strictly hyperbolic. Denote 



by Rpf[g], A]y[g} the eigenvector and eigenvalues of Dg that are associated with the N- 
characteristic family, and define 7rj^ [g] (u) , . . . , 7rjy [g] (u) in the usual way. To exhibit a 
state satisfying (12.441) we consider the following vector field W[g] : Bs - 



9 N. 



( [#] («!,■■■ 



W[g]( Ul 



UN 



\ 
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From (I2.41j) - (j2.43|) one sees that, for any (ui, . . . , u^) G Bs, 



N-l 



W[f]{ Ul , ...,u N )- K . . . , u N ) = + !) ! «? + (iV + 2)! < + P( Wl , . . . , 



mat), 



where P(ui, . . . , Wjv) is a polynomial expression in iti, . . . , un such that which each term is 
of degree 3, at least. It follows that for some 5 > and for all 5 G (0, S ) and (ui, . . . , u N ) G 
Ss (the sphere with radius 5) 



In consequence, for 5 G (0,£o), the index of the vector field W[f] around the sphere Ss is 
1 (by a standard homotopy argument, for instance). Now, it follows that for all g close to 
/, the vector field W[g] does not vanish on Ss and is of index 1 around S$. Hence, W[g] 
necessarily vanishes inside Bs, and this establishes the claim. □ 

3 A second proof of the inner speed variation esti- 
mates 

3.1 Notation and preliminaries 

In the present section, we rely on an alternative approach to the Riemann problem due 
to Bianchini and Bressan [HI El E] and we provide a second proof of the inner speed 
variation estimate already derived in the previous section. 

Throughout, we consider a strictly hyperbolic system of conservation laws ( 12. lft and 
use the notation of Section [2j Note that the parametrization (often denoted by the letter 
s) along the wave curves will be here different from the one in the previous section. The 
sup norm in M. N is used throughout the present section. 

The following framework is based on a prescribed family of "traveling waves" ; for the 
motivations of which we refer the reader to (23, UHl EH [8], E] and to the discussion in 
Section H] below. Consider a given family of smooth, vector- valued maps rj = rj(u,Vj,aj) 
for (u, Vj, <Tj) G R N xRxl (presumably associated with viscous profiles of a regularized 
version of (12.11) ). Given a base point which can be assumed to be the origin in Mr, and 
setting $ := lj(0), etc., we can normalize the vectors rj so that 1°, ■ rj = 1. 

From rj we can determine the speed functions Xj 



...,u N )- (u u ...,u N )>0. 



Xjfav^ffj) := lj ■ 



and we assume as in [H] that, for some constant Co > 




(u,Vj,aj)\ < C \v. 



(3.1) 



do-j 
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and 



d\ ■ d\ ■ 

3 (u,Vj,aj)\ <C \u\, \-^-(u,Vj,aj)\<C Q \u\\vj\. (3.2) 



dvj 1 ' 1 daj 



In particular, it follows that \j(u, 0, Oj) = Xj(u). 

Then, the wave curves associated with the system (12.11) are constructed in the following 
way. Fix some small Si < 5o and, for j = 1,...,N, s G [0, 6i), and u_ G B$ , define a 
family of curves Tj(s; u_), issuing from u_ and of "length" s, as follows: 



r,( S ; W _) := | 7 (r) = (uW.^W^^t)) G LiKMjR** 2 )/ 
m(0) = U-, fij(u(r)) = [ij{uJ) + r, 
^•(0) = 0, \ Vj (t)\ < 8 U |a,-(r) - A°| < 2C ^ 

where := 1® ■ u. To any curve 7 G Tj(s; u_) we associate its j-flux function 

fjh}( T ) ■= / ^(u^ffjOCO^- 



Define also the nonlinear operator Tj : 7 1— > (u,Vj, dj) on rj(s,M_) by 

u(r) = u_ + / rj(u,Vj,aj)(r') dr' , 



= /jHtj) - c o nv [o,.]/j[7](^), ( 3 - 3 ) 
^iO") = ^( conv [o,s].S[7](^))- 
Define a distance between 7,7' G rj(s,M_) by 

-0(7,7') : = - m'IU°° + Ikj - UylU 1 + ll^-cTj - v 'j°'j\W- (3-4) 

It was established in [7] (Proposition 3.2 therein) that, for b\ suitably small, the operator 
Tj is a contraction with constant 1/2 in the metric space Tj(s;u-) endowed with the 
distance D. To any s and u_ we can thus associate a unique curve 7* = r y* su _ that satisfies 

This allows us to introduce the wave curve ipj = ipj(s; tt_) by keeping, for every s > 0, 
the end-point of the curve at r = s, only, i.e. 

il>j(s]uJ) :=7,V( S )- 

For every s > 0, the state w_ is connected to the right-hand state u + := ipj(s; uJ) along the 
curve 7* = (w*, w*, (Tj) (associated with the given s, uJ) and we can define the corresponding 
j-wave fan: 

x/t < (Tj(0), 

u(t), x/t = <Jj(r), 

U_, x/t > <Tj(s). 
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When s < a similar construction is done using the concave hull (denoted below by cone) 
instead of the convex hull. 

The wave curves are Lipschitz continuous and, by the implicit function theorem, any 
Riemann problem (ui,u r ) can be solved uniquely, by combining wave curves together, 
i.e. 

u r = \P(s; ui) := ^ N (s N ; ifi N -i(s N -i, ■ ■ ■ «i)))> 

where we use the notation s := (s x , . . . , sjv)- 

To estimate the interaction between waves we introduce Bianchini-Bressan's amount of 
interaction [7j as follows. First of all, consider two waves of the same family 



u + = ^j(s,u-), u'+ = i/)j(s';u'_), 



(3.5) 



together with the corresponding curves 7 and 7' and the corresponding j-flux fj = fj[j] 
and fj = fj[j']- Define the real X = X(it_, u + ; u'_, u' + ), as follows, where for definiteness 
we assume that s > 0. (When s < one should replace all convex/concave hulls by 
concave/convex hulls.) 



If s, s' > 0, we set 



rE[0,s], 



fj(s) + f'(r-s), tG[s,s + s'}, 



and 



conv [0iS ](^ U - conv [0iS+s1 (/J- U 



+ 



conv 



[8 , a+ M u /;)(o - conv^+y]^- u /;)(o 



If — s < s' < 0, we set 

X(u_, u + ; it'_, u' + ) 



s+s' 



s+s' 



conv[ 0| ,]/j(0 - conv[ 0|S+ ^/i(0 
conv[ , s ]^(0 - conc [s+s /, s ]/i(0 



(3.6) 



If s' < — s < 0, we set 



conc [s i,o]/j(0 - conc [s / ,- s ]/j(0 



concr,/ 



.0]/j(0 



conv [ _ Si0] /j(0 
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More generally, define the amount of interaction associated with two Riemann problems 
u r = tf(s;u,), u' r = *(s';«0 by 

T{uuU T \v! h v! r ) := ^2 \ s i s 'j\ + ^ 2 "( tt -) tt +; tt -) tt +) i 

where the second sum is over all j-waves (u-,u+) and (u'_,u' + ) in the Riemann problems 
(ui,u m ) and (u m ,u r ), respectively Now the interaction functional Q is defined as 

r s 3 r s 'j 

Q{ui,u r ,u' l ,u' r ) = ^2\s i s' j \+22 / \aj(T) - trj(r)|. 

It will be convenient also to set Q(/Uj, it m , w r ) := Q(ui, u m \ u m , u r ) for all ui, u m , u r . 

Finally we will need another distance between curves. Consider two curves 7 G Tj(sj \ u) 
and 7' G R^s^w'), and restrict attention to the case where Sj and s'j have the same sign. 
If, for instance, both of them are positive we set §_j := min(sj, s'j) and 

P(l,l') ■= Dil^hl'io,^) + \ s j ~ s'jl 

where the distance D introduced earlier is extended in an obvious way to any two curves 
with the same length but not the same base point. 



3.2 Interactions estimates on the inner speed variation 

Recall from [7] the following version of Glimm's interaction estimates. 

Theorem 3.1. 1. Consider any three states Ui,u m ,u r such that u m = ty(si m ;ui), u r = 
ty(s mr ;u m ), and u r = ty(si r ;ui). Then, the corresponding wave strength vectors satisfy 

7TLV I 

l(ui (3.7) 

2. Let 7j m , 7™", and 7j r be the curves associated with the j -waves in the Riemann 
solutions under consideration. Define jj by 

7™J J (• - 4 m ) on [4 m + sf r > 0], s™ r < -sf < 0, 
where 7^ m U 7™" zs defined by 




Ira ,j w . = J fj v < r < sf, 

'1 'J I ^,lrn( „lm\ _i_ mr/ e '»rA c im <^ t <^ e' m -I- e mr 

~r / .7 v & j J' & j — — j i ' 
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(the definition when s l J n ,sJ lT are negative being similar). Then, the curves satisfy the 
following interaction estimates: 

3 

3. Consider a piecewise constant function with small total variation u : IR — > R . Let u' 
be the function obtained from u by replacing two consecutive Riemann problems (uj.,«(. +1 ) 
and (Mfc + i,Mfc +2 ) by Uk+2)- Then, for some c> we have the estimate 

Q{v!) < Q{u) - cl(u k , u k+2 ). (3.9) 

4- In the situation above, let Q' be the interaction amount associated with u but with 
the i-th curve Yk ^ n the Riemann problem (uk,Uk+i) being replaced by 7^ corresponding to 
another problem u[ = ijji(u' ). This modifies the Q in the following way: 

\Q'-Q{u)\<P{ihi^. (3.10) 

We now turn to the investigation of the inner speed variation (defined in the previous 
section) which here takes the following form: given u r = iftj(s, u{) with corresponding curve 

"&j{ui\ U r ) = O-j(s) - CTj(O). 

We provide a new proof of the estimate already derived in the previous section. 
Proposition 3.1. (Inner speed variation - same family.) Consider two j-wave fans 

and denote by 7; m = (u lm , vl™ 1 , al" 1 ) and j mr = (u mr ,v™ ir ,aY ir ) the associated curves. Set 
u 'r := ^j{ s im + s mr'i u i) an d denote by 7' := (u',v'j,o~j) the associated curve. Then, the 
inner speed variation of the j-wave fan (ui,u r ) satisfies the following estimates (assuming 
Sim > 0): 

• Monotone case : < si m ,s mr . 

&j(ui,u' r ) <max (i9 j (ui,u m ),^j(u m ,u r )) 

+ (of (0) - of (0)) + + 0(1) X{u h u m , u r ). 

• Non-monotone case (I) : < si m + s mr < si m . 

M u i><) < ^j(u h u m ) + 0(l) \s mr \. (3.12) 

• Non-monotone case (II) : sj m + s mr < < si m . 

r d j (u h u' r ) < tijiu^Ur) + 0(l)s lm . (3.13) 
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Moreover, completely similar estimates hold in the case s/ m < 0. 

Proposition 3.2. (Inner speed variation - artificial wave fronts.) Consider three states 
Ui,u m ,u r satisfying solely 

and let 7 = (u,Vj,aj) be the curve associated with the wave (u m ,u r ). Introduce u' m : = 
ipj(s;ui) and the corresponding curve 7' := (V, v'j, a'j). Then, it holds 

&j(ui,u' m ) = fij(u m ,u r ) + 0(s) \u m - Ui\. 

A completely similar statement holds when (ui,u m ) is a j-wave and (u m ,u r ) is an arbitrary 
jump discontinuity. 

In fact, from the above statement, we can also deduce the desired estimate for the 
interaction of waves of different families. 

Proposition 3.3. (Inner speed variation - different families.) Consider three states ui, u m , u r 
such that, for i 7^ j , 

Urn 1pj\Slm,i Ul) , U r Ipii^Smn " m ), 

and denote by 7'™ = (u lm , v 1 ™, a 1 ™) and 7™" = (u mr , w 4 mr , a™ r ) the curves associated 
with the wave (ui,u m ) and (u m ,u r ), respectively. Set also u m := xfji(s mr ;ui) and u r := 
^i{si m \u m ), and denote by i lm := (u' lm , v'\ m , a'] m ), i mr := (u' mr , ^' mr , a'J nr ) the corre- 
sponding curves. Then, one has 

$j{Um, u r ) = i}j{ui, u m ) + 0(s lm s mr ), 

$j{u h U m ) = 1?i(li m , U r ) + 0(si m S mr ). 

Corollary 3.1. Consider four states Uo,u' ,Ui,u[ such that 

ui = if)j{s\ M ), u[ = ipj(s'; u' ) 

with < s < s' . Denote by 7 = (u,Vj,o~j) and 7' = (u', v'j, er'-) the curves associated with 
the waves (u ,ui) and (u' ,u[), respectively. Then, one has 

- oj(0) - aj(s) + ^(0)1 < \u - u' \ + |a' - s\. (3.14) 

Clearly, Proposition 13.31 follows from the fact that 

s < \^j{s;u) -u\< S. 

We postpone the proof of the other statements to Subsection 13.41 In the next subsection, 
we derive some additional properties of the Riemann solver constructed via (13.31) . which 
will be useful in our proofs but are also of independent interest. 
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3.3 Additional properties of the Riemann solver 

To derive the inner speed variation estimates stated earlier, it will be convenient to de- 
compose the problem of finding a fixed point of (13.31) in two steps: on one hand, finding 
the geometric wave curves (that is the component it) and, on the other hand, finding the 
speed of propagation of the wave fan (that is the components (vj, o~j)). This is the subject 
of the following proposition. Consider the set 



T,( S ;«_) := {(O.o-jO) G Lip([Q, s]; M 2 ) / 

Vj (0) = 0, \ Vj (r)\ < S u \aj(r) - A°| < 2C 5 1 

endowed with the sup norm \\(vj, fj)||oo : = Halloo + Halloo- Given u* G Lip([0, s];R N ) 
satisfying u*{r) = u_ + r and \u*{t) —U-\ < r, consider the operator 



a at 



3\ J ~ dr 



conv[ 0)S ]/j[tt*, vj, aj] (r). 



Lemma 3.1. Fix u* G Lip([0, s]; M. N ) satisfying u*(r) = u_ + t and \u*(t) — u_| < r. 
Provided the range of u* is included in a sufficiently small neighborhood of 0, the operator 
Q" admits a unique fixed point in Tj(s; uS), which will be denoted by {V^ , E" ) . Moreover, 
for some C\ > 0, 

\\(V s u \^f) ~ (V s u \^f)\\oo < Ci \\u l -u 2 



,2|| 

oo- 



Proof. We will prove that if the range of u* is included in a sufficiently small neighborhood 
of 0, then the map is a contraction with constant 1/2. Indeed, one has 

||fif (Vi,<7i) - ft** (V 2 ,02)||oo ~ ll^/j-("j u *J U l> CT l) - ^rX(->W*>«2,ff2)||oo 

< ||Ay(u*,Ui,cri) - (li* , , cr 2 ) 1 1 oo . 

With (JS2D this yields 
||ftf (v 1 ,a 1 )-nf(v 2 ,o- 2 )\\ oc 



< 



C \\u* - tt ||oo(l + IK'||oc)|IOi,o"i) - (v 2 ,a 2 ) 



This proves the first part of the lemma. 

To establish the second part, we use the above contraction property and write 

ikk^e^-o^s^iu < ii^(K U2 >sr)-(K W2 ^r)iu 

< IK-U2II00, 

where the regularity of the function A,- has been used. This completes the proof. □ 
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We now state some properties of the curve describing the wave fan, under the as- 
sumption that the w-component is already known. Introduce the translation operator 
t s : / h-> T a f, defined by 

r s (f)(y):=f(y-s). 

When a function u* is defined on an interval larger than [0, s], in order to simplify the 
notation, we simply write (V,"*,!!"*) instead of (V s E S I[M ). 



Lemma 3.2. (A splitting property.) Let < si < s 2 and u* G Lip([0, s 2 ], K"^) 6e as m 
Lemma \3. 11 7/ 

Vtf(ai) = 0, 

iaen (V^, 2 *,E"*) coincides with (V^*,E"*) on £/ie interval [0, Si], and coincides with the 
translate ^(V^-^ , ^Ij-si ) on ^ e ^nierva/ [si,s^. 

Proof. It suffices to check that (V^ , E"J)|[ 0;Sl ] and r_ Sl (V^*, ^s 2 )|[si,s 2 ] are fixed points of 

the operators Qs[ ' and f2 S2 _ Sl , respectively. But this property is a direct conse- 
quence of the following fact: for every Lipschitz continuous function g : [0, s 2 ] — > M and 

si e [0, s 2 ] 



,conv 



[o, S2 ]#)Oi) = #00 conv [0iS2] o 



conv [0jSl ]O on [0, si], 
conv [si)S2 ]5 on [si, s 2 ]. 



□ □ 

Now, we prove (compare with P, Lemma 3.3]): 

Lemma 3.3. (A monotonicity property.) For any < si < s 2 and -a* G Lip([0, s 2 ], M^) 
as m Lemma l3J[ it holds: 

E£>E£ on[0, Sl ]. (3.15) 

Proof. Set 

s D :=max|sG [0, Sl ], K 2 * (s) = o}, 

and note that the function conv[ 0iS2 ]/j(r, u*, V^* , E"*) is affine on the interval [s D , si]. Sup- 
pose first that s D > 0. We have V S2 *(s D ) = and, by Lemma l3~2l the curves (V S "*,E"*) 
and (V S2 *, E"*) coincide on the interval [0, s n ]. 

Hence, it is clear that (13.1 5ft is valid in [0, s D ], so that on the interval [s n , si] we obtain 

s:;( 5 )>s:;( S d ) = e:;(s d ) = e:;(s). 

Recall here that, since 7 is Lipschitz continuous, the function fj^) is W 2,oc . 

Suppose next that s n = 0. Because conv[ 0)S2 ]/j(r, u*, V™* , E" 2 *) is affine in [0, si], it 
suffices to establish (13.151) at s = 0. We distinguish between two cases: 
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either is an accumulation point of values s such that V™*(s) = 0. In this case, we 
have 

and, on the other hand, 

e£(o) < Ai(u*(0)) = ^(0,«M^*,EJ). 

or s A is the smallest positive r such that V^*(r) = 0. In that case, both £"* and 
£"* are constant in the interval [0, s A ], the former coinciding with in [0, s A ] 
thanks to Lemma 13.21 Now, let us use (13. 2p and recall that both and E" 2 * 
are constant in [0, s A ]. Recall also that both V u & and V u 2 are linear in [0, s A ] and 
V s u 2(s A ) = V£(0) = V£ (0) = 0. We then find: " 

\\%{u\ v;:,^l) - %{u\ V£ , K:)\\l~(o,s±) 

< \\Xj(u*, V^*,£"l) - Xj(u*, V^jf , E"*)|| L i(o )S A) 

< \\u*\U(V s t,^l) - (V£,K:)\\l H o, s ±) < 

at least if u* remains within a small neighborhood of the base point 0. This yields us 



which leads to (13.151) and completes the proof of Lemma 13.31 
From Lemma [3.31 we deduce: 



□ 



Lemma 3.4. Let s 1 , s 2 > 0, s e (0, si), and w* G Lzp([0, si + s 2 ]; M^). // 

VS;„(S) = 0, (3.16) 

then 

Kf(i) = 0. (3.17) 
Similarly, if (13.161) holds for some s G (si, Si + s 2 ), i/ien V^ -81 " (s — si) = 0. 
Proof. We only prove the first statement since the second one is similar. We introduce 

U -~ T - sU \[s,s 1 +s 2 ]i V s 1 +s 2 -si Z ^s 1 +s 2 -s- 

Using Lemma [3.21 we obtain 

= — convp^+^/iCs, «*, Sj +Sa , Kf+s 2 )- 
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Using V™* +s (s) = 0, we then have 

= A i («*(S),0 ) E?; + „(S)) = A i (u*(S))- 
Now, from Lemma 13.31 and 13.21 we deduce 

r ? E*_~(3) > r,S: i+S2 _,(?) = S«* +S2 (5), 

It follows that 

i*K-&) > A!(«*(S)). 

But, on the other hand, 

r^_,(i) = |-conv [ o, Sl _ 3] /,(«,K , :- ? ,S: i _ ? )(0) 

= A 1 («*(S),0 ) Ej_ y ) = A 1 («*(i)). 

Hence, we obtain 

7 ^ E W"i(S) = A 1 (tt*(S))=^(S). 

But, we have also 

^1^-^(5) = o = v£^(g), and 7sfi||p^_s,(S) = u*(3). 
We introduce the curve (V, E) : [0, si] — > M 2 given by 

MeK^. ^ (3.18) 

Clearly, 7 is Lipschitz continuous. To prove (13.171) . it suffices to prove that (V,"*, = 
(V, E), and hence that 7 is a fixed point for fi" . 

Now the following property is satisfied by any Lipschitz continuous function g : [0, si] — > 
JR.: if its convex hulls /ii := convro,ii<7 and := convpisilfl 1 satisfy 

/ii(s) = h 2 (s), h[(s) = h' 2 (s), 

then conv[ 0iSl ](7 coincides with /ii in [0,3] and with /12 in [3, si]. Using the definitions of 
(V^*,Ef ) and (V£_ y ,E*_ y ), one deduces (V™* ,EfJ = (F,E) and hence 1ETIB]) . This 
completes the proof of Lemma 13.41 □ 
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We conclude this section with a technical observation which will be useful in the course 
of the proofs of the above statements. 

Lemma 3.5. All of the curves 7 obtained by a fixed point argument with the operator (13. 3p 
are Lipschitz continuous with a uniform Lipschitz constant. 

Proof. Consider first the component (u, Vj) part; the desired estimate is elementary in 
view of the property ^(conv/)^ < ||<i/||oo- It remains to consider the <x,-part and, 
precisely, to prove 

MT)-<T j (/)\<C\T-7'\. 

First of all, observe that for every / £ W 2 '°° and < x < y < s the convex hull 
h := conv[ 0s ]/ satisfies 

osc(h') < osc(f'). 

[x,y] [x,y] 

(This is easily check by considering the contact points between / and h in [x, y]). Hence, 
we have 

and it follows that 

osc(o-j) < II^AjHoo osc{u) + IISuAjlloo osc{vj) + H^AjHoo osc(o-j). 

[t,t'] [ t : t '] l T : T '] [ r i T '] 

Now, we already have uniform Lipschitz estimates on (u, Vj), and it follows from (13.21) that 
the coefficient in front of osc(aj) is less than 1/2, which allows us to conclude. □ 

3.4 Proof of the inner speed variation estimates 

We begin with an elementary lemma. 

Lemma 3.6. Given g £ W 2 '°°([—M, M]) , there exists a positive constant C\ = C\{g) 
(depending upon M and sup \d 2 g/du 2 \) such that, for all —M < a < b < c < M, one has 

\\(conv [a>c] g)\ [ab] - (conifaqgyWoo < G x (c- b), (3.19) 

\\{conv [ajC] g)\ [bc] - (conifoc^yiloc < C t (b - a). (3.20) 
A similar statement stands for the concave hull. 

Proof. We limit our attention to the proof of (13 . 1 9[) in the case of the convex hull since 
the other inequalities are proven in the same way. We use the notation h\ a u := convr a M<7, 
etc. In view of this definition we note that 

h[ a ,c] < h[ a>b] < g in [a, 6], 
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hence we can restrict attention to the case where a and c are the only contact points 
between g and h[ a>c ]. Hence one gets 



, _ g{c) - g{a) 
^ = c-a ■ 

^ s h'\ a c ] (a) < h'< a 6 i (a) it is sufficient to prove that 

. / /.n ff(c) - g(o) 

^m(^) — — ; — < 0(c - b). 

Suppose that h\ a ^ and g are not tangent at b, and denote by d < b the "last" contact 
point of h\ a u with g. It follows that 

>/ /n g(g) ~ g(g) _ g(d)-g(b) g(a) - g(c) 

h wiW- c _ a - d _ b - 



< 



a — c 

^(d) - p(6) ^(rf) - g(c) 



d — b d 



In the latter inequality we used the fact that the point (d, g(d)) is above the line connecting 
(a,g(a)) to (c,g(c)). The latter term in the inequality above is obviously O(o — b). 
Now, in the case where h\ a fi, and / are tangent at b, the difference is 

m _ 9(c) -9(a) = gl{h) _ 9(a) -9(c) 



< 9'(b) 



a — c 
9(b) -9(c) 



where we have used that (b,g(b)) is above the line connecting (a,g(a)) to (c,g(c)). The 
right-hand side is also clearly 0(c — b) and this completes the proof of (13.1 9ft . □ 

Proof of Proposition 13.11 To simplify the notation, we use the index 1 (resp. 2) to 
denote the objects (s, 7 = (u, v iy <Tj)) relative to the left-middle pattern (that is, Im) (resp. 
the middle-right pattern mr). We drop the index i corresponding to the wave family 
(which is the same for all objects considered here). Recall that a prime indicates an object 
corresponding the outgoing 2-th wave. 

We distinguish between three cases: Si,s 2 > 0, —s\ < s 2 < and s 2 < —s\ < 0. The 
cases where si < are treated similarly using the concave hull. 

Case 1 : Assume s\, s 2 > 0. 

Let us define 71 U 72 by ( 13. 8ft . let be the u-part of it. It follows from Theorem 13.11 and 
Corollary 13.11 that 

D(y, 7! U 72) = O(l) X(u h u m , u m , u r ), 
where the distance D is given by (13.41) . In particular, 

|| M # - M '|| 00 = 0(1)X, 
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(3-21) 



which in turn yields, using Proposition 13.11 

\\(v s f +S2 ,Ef + j - (Kf +S2 ,s:; +S2 )iu = am. 

Hence, to establish (13.111) it is sufficient to prove 

< max (a\ Sl ) - a\0), a\s 2 ) - a 2 (0)) + (a 2 (0) - a 1 (0)) + . 
Let us introduce 

s a : = max js G [0, si], K"t S2 ( s ) = °}> 

s 6 := min js G [s l7 s x + s 2 ], K"ts 2 ( s ) = °}- 

From Lemmas 13.21 and 13.41 we deduce that H™* +S2 coincides with in [0,s a ], and with 
t^E"! in [sb, si + S2]. Moreover, is constant in [s a , s&]. Therefore, we see that 

• if s a > and Sb < s\ + S2, then 

S:f +S2 ( Sl + S2 )-E:f +S2 (0) = £-( S2 )-££(0) 

= e: 2 2 ( S2 )-e: 2 2 (o) + s: 2 2 (o)-s^(o), 

which yields (I3.2ip . 

• if s a > and s?, = si + s 2 (the equivalent could be done in the case s a = and 
s b < si + s 2 ), then £ S1 + S2 is constant in [s a , Si + s 2 ], hence 

Srr +S2 ( Sl + S2 )-Sf +S2 (0) = E*(a 8 )-^(0), 
which yields again f)3.2ip . 

• if s a = and Sb = si + s 2 then S" 1+S2 is constant in [0, s\ + s 2 ], hence the result is 
satisfied. 

Case 2 : Assume — si < s 2 < 0. 
Let us write: 

1* ■= 7l|[0,«i+« a ]- 
It follows from Theorem 13.11 and Corollary 13.11 that 

hence 

\\(vf +S2 ^f +S2 ) - (Kt+^+jiu = o(i)i. 
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We note that, as all the speeds are bounded, X = 0(s 2 ). Hence, to establish (13.121) it 
suffices to show 

cr*( Sl + s 2 ) - a*(0) < <n(si) - oi(0) + 0(s 2 ). 
As previously, using the contraction property, it is sufficient to prove that 

\\^* + s 2 ( v U a l)\[0, Sl +s 2 ] ~ (Vl,0l)|[O, Sl + S2 ]||oo = 0(S 2 ). 

Computing the difference, we are led to establish that 

||[conv [0iSl] /(7 1 )]| [0iSl+S2] - conv| [0iSl+S2] /(7 1 )|| £ i P = 0(s 2 ). 

Using Lemma l3~ol the expression of fi and the smoothness of Aj, we see that it is sufficient 
to prove a uniform Lipschitz continuous estimate for the curves 7. This is precisely given 
by Lemma 13.51 



Case 3 : Assume s 2 < — si < 0. 

This case is done similarly as the previous one, and this completes the proof of Proposi- 
tion O □ 

Proof of Proposition 13.21 In the sequel, all objects (s, 7 = (u, Vj, <7j)) without a prime 
are relative to the incoming j'-wave, the objects with a prime to the outgoing one. Without 
loss of generality, we assume that s > and use convex hulls. We treat only the case that 
the given j-wave is (u m ,u r ), since the other case is similar. Define 

u(-) :=«(•) -u m + u h 7(-) := (u, v, ct)(-), 

and consider d = u m — u\. For < r < s we have 

~ d ~ d 

Xj(j)(s) - A^xKr) = —Xj(u(s)+<;j^,v(s),a(s))d<;. 

In consequence, using the uniform bound available on the second-order derivatives of Xj, 
we find 

l|Ai(7)(0 - Ay(7)(0 + {Ai(7)(0) - A i ( 7 )(0)}|| i co (M) = 0(s\d\), 
and, in consequence, 

\\m)(r) - 7,(7)(r) -r{A,(7)(0) - A,( 7 )(0)}||^ = 0(s\d\). (3.22) 

Next, since ||conv/ — convg||vyi,oo < ||/ — we deduce that 

llconv^]/^) -conv [0l .]^(7) -r{Aj(7)(0) - A^XO)}^,- = 0(s\d\). (3.23) 
Let us now show that 

||(V* E*) - (v,a + {A,(7(0)) - A i ( 7 (0))})|| oo = 0(1) s\d\. (3.24) 
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In view of (13 .ip we have 

A,(7(0)) = A J (n m ), A J -(7(0))=A i (« I ). 

We denote a := a + A 3 -(7(0)) — Aj(7j(0)) and 7 = (u, v, a). Using the contraction property, 
to establish ( 13. 24ft it suffices to show 

\\(v,a)-to%v,a)\\ 00 = 0(l)s\d\. 

Recall that ||dconv/ — dconvgW^ < \\df — dgWco, and observe that g — convg remains 
unchanged if one adds an affine function to g. Recalling also the definition of the operator 
fi", we obtain 

= Wfjil) - conv [0ja ]^(7) - J)(l) + conv[o, s ]£-(7)||oo 
+ II^I conv M/j(7) + Aj(7(0)) - Xj(jj(0)) - —conv [0^^(7)1100 

£ II fj (7) - conv [0iS ]/ i (7) - fj(j) + convp^/y^Hoo 
+ \\fj(l) ~ conv [o, s ]/i(7) - fAl) + c o nv [o, s ]/j(7)l|oo 
+ 11^(7(0)) - ^(7(0)) + ^/i(r,7) - ^/;(t,7)Hoo 

+ II ^/i( r >7) - ^:/i(r,7)||oo. 
Taking ( I3.22l) -( l3~23l) into account, this yields 

||0,5-) — ^"(V, Cr)||oo ^ \\\j{T,U,V,&) ~ \j(T,U,V, d-)\\oo + 0(1) s\d\, 

and, using (13.21) . 

||(«,cr) < ||t;||oc|k - c^Hoc + O(l) 

< 0(i)|MUw m -^| + o(i) s |d|. 

Since i> = 0(1) s (thanks to the bound on Xj in the domain under consideration) we 
deduce (]3.24|) . and it follows that 

£ s (s) - £"(0) = S u (s) - S u (0) + 0(1) s|d|. 

In view of the proof in Bianchini [7] (cf. Lemma 3.8 therein) we have 

D{z i ,i)=0{s)\u m -u l \. 

(This follows by the same procedure, it is sufficient to prove that .0(7, T l (7)) is of order 
0(1) |s||« m — Ui\.) From Proposition 13. II we deduce 

\\(v u \j: u ')-(v^^)\\ 00 = o(i)s\d\. 
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which yields the desired conclusion. This completes the proof of Proposition 13.21 □ 

Proof of Corollary 13.11 In view of Proposition 13.21 we see that it suffices to treat the 
case u = u' . One can estimate Hu 1 — ^r J|l°° by 

h ~ u '\[o,s]h°° ^ D (V(l\[0,s]),l'\[0,s])- 

Recalling Lemma [3.61 and the uniform W 2 '°° estimates on fii'j 2 ) (by the regularity of Aj 
and Lemma l3~5|) . we deduce that 

w ([o, s ]IU°° ~ I s ' - 4 

In turn, with Proposition 13. 1[ this yields — E^'H^oc < \s' — s\. Using the uniform 
Lipschitz continuity estimates on <jj, to prove (13.141) it suffices to check 

||(E«')| [0 , S] -E«'||^<| S '- S |, 

which, again, follows from the contraction property and Lemma 13.61 □ 

4 A generalization to nonconservative hyperbolic sys- 
tems 

We now turn our attention to nonlinear hyperbolic systems in nonconservative form, i.e. 

d t u + A(u) d x u = 0, u = u(t,x) e R N , t > 0, x e R. (4.1) 

As usual, we assume that the matrix A(u) has real and distinct eigenvalues Xj(u) and a 
basis of left- and right-eigenvectors lj(u),rj(u), normalized as in Section [21 Our aim is to 
extend to the nonconservative system (14.11) the theory of the Riemann problem discussed 
in previous sections. Recall that the distributional definition of solution does not make 
sense for nonconservative systems, and that a suitable notion of weak solution for such 
systems was introduced in [2EI EH EI] which is based on a prescribed family of Lipschitz 
continuous paths. The Riemann problem for genuinely nonlinear and nonconservative 
systems was solved in [26], [16] . We are here interested here in the generalization to systems 
that need not be genuinely nonlinear. 

Following LeFloch [271 El] we can define the family of paths and, therefore, the gen- 
eralized Rankine-Hugoniot relation for (14.11) from the family of traveling wave solutions 
associated with a regularization with small diffusion 

d t u + A(u)d x u = ed x (B(u)d x u), (4.2) 

where the diffusion matrix B = B(u) is, say, positive-definite. 

We first observe that Bianchini-Bressan's arguments generalize straightforwardly to 
nonconservative systems. More precisely, the Riemann problem admits a unique entropy 
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solution if the maps rj satisfy the conditions imposed in the beginning of Section [SJ On 
the other hand, the interactions estimates (stated in Theorem 13.11) were established in [H] 
in the case that the nonconservative system is regularized by an identity diffusion matrix 
B(u) = Id. The rest of this section will thus be focused on extending Iguchi-LeFloch's 
method. 

Using the notation in Section [3] introduced within Bianchini-Bressan's method, this 
amounts to fix vectors Tj and to introduce the j-Rankine— Hugoniot curve Hj(s; uJ) = 
u(s,s;u_) together with the shock speed Xj(s, it_) = <jj(s, s;u_), where we define the 
map r \—> (u, Vj, Cj)(r, s; uJ) by the differential system: 

d T u(r) =r j (u,v j ,a j )(r), 

d T Vj(T) = Xj(u, vj, aj)(r) - <r 3 -(0), , 4 ^ 

1 f'~ 

(Tj (t ) = crj(O) = - \j(u,Vj,(jj)(t)dt, 
s Jo 

with initial condition 

n(0,s) = n_, Vj(0, s) = 0. (4.4) 

Theorem 4.1. When the shock curves are defined by (14. 3p . all of the results in Section^ 
i.e. the Riemann problem, regularity of the wave curves, and interaction estimates for wave 
strengths and inner speed variation, remain valid for the nonconservative system 04. ip . In 
particular, the Hugoniot curves defined by (14.31) satisfy 

sd s \j{s]uJ) = Kj(s;u-) (\j(Hj(s;u^)) -\j(s;u-)), 

d s Hj(s,u-) = rj(Hj(s,u-)) +Kj(s;u-)(Xj(Hj(s;u-)) - Xj(s;u-)), 

where the function Kj = Kj(s\uJ) > is smooth, bounded, and bounded away from zero, 
and satisfies Kj{~s; uJ) = 1, and Kj is a smooth vector-valued map. 

Observe that no assumption (of genuinely nonlinearity) need be imposed on the ma- 
trix A. Our result extends the construction of the Riemann problem given by Dal Maso, 
LeFloch and Murat [H] for genuinely nonlinear systems. Dealing with non-genuinely non- 
linear systems is more involved. 

Observe also that no reference to the vectors Tj(u, Vj, o~j) is made in (14.51) . so that some 
further generalization of the above theorem is possible; see Remark 14.11 below. 

Proof. The analysis in Section [2] relied on the key property (12. 3p of the shock speed and 
characteristic speed along the Rankine- Hugoniot curve. It is remarkable that this property 
remains true when the shock waves for nonconservative systems are defined from traveling 
waves. Differentiating (14. 3p with respect to s, we obtain a linear differential system for the 
"unknowns" u' := d s u, v' := d s v, and a' := d s o~: 

d T u' = V u rj(u, Vj, o~j) ■ u + d v Tj{u, Vj, Cj) v' + d a rj{u, Vj, o~j) a', (4.6a) 
d T v' = V u Xj(u, Vj, o~j) ■ v! + d v Xj(u, Vj, o~j) v' + d a Xj(u, Vj, o~j) a' — a', (4.6b) 
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1 1 f s r ~ 

°' = -(^j( u ( s i s ~, u -)) - a ( s )) + ~ \ V u \j(u, Vj, (Tj) ■ u'(t, s) 



+ d v \j(u, Vj, (jj) v'(t, s) + d a \j{u, Vj, aj)a'(s)^dt. (4.6c 



The term Xj(u(s, s; u_)) in f)4.6cl) arises from the fact that v j(s, s, uJ) = and \j(u, 0, ■) = 

In view of g3D, we have u'(0,s) = and v'(0,s) = 0. Hence, from (Oall - ljPKjl it 
follows that, given any s, 



'|U°°(0,r) + |K ;/ ||l°°(0,t) < 7-(||w'|U°°(0,r) + L°°(0,r) + (4.7) 



Hence, using ( 14.71) within (j4.6cj) and noting that, by fl3.ip - fl3.2l) . we have \d a Xj(u, Vj, Oj) \ « 
1, we can deduce (for sufficiently small r, s) 



\d s a(s) 



< 



-(Xj(u(s, s; U-),Vj(s, s; «_), cr(s; u_)) - a(s)) 



(4.8) 



In other words, we have established the first equation in f)4.5p with Kj bounded. 
Considering next the map Hj we see that 

d s Hj(s;u-) = d T u(s,s;ii-) + d s u(s,s;u-) 
= fj(u,v,a)(s) + dgU(s,s;uJ). 

To handle the first term, we observe that Vj(s, s; it_) = and fj(u, 0; a) = Tj(u). For the 
second one, we use again f)4.7p which, together with f)4.8p . yields the second equation in 
(14. 5 p with Kj bounded. 

We next establish some regularity on Kj,Kj- First, from (14.31) it follows that, for every 
fixed s, the mapping (u,v,a) is of class C°° with respect to the variable r. On the other 
hand, it was proven in [H] (see (3.26) therein) that this mapping is Lipschitz continuous 
with respect to the variable s. Now, consider fl4.6al) - (14.6bp as a linear differential system, 



with Lipschitz continuous coefficients, the terms cointaining d s a being viewed as a source- 
term. (Recall, moreover, that d s a is independent of r.) Then, it follows that v! and v' can 
be written in the form 

(«', v')(r, s) = d s a(s) (h u (r, s), h\r, s)), (4.9) 

where h u , h v are Lipschitz continuous in both variables and 0{r) at most in the sup norm. 
Taking (14.91) into account into (I4.6cl) yields 

/ \j{u(s,s\u-)) - a(s) 
°~ s s(l-p) 

P ■= - (VuXj(u,Vj,aj) ■ h u (t,s) 



s 







+ d v \j(u, Vj, fjj) h v (t, s) + d a \j{u, Vj, (Tj) ) dt 
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where the coefficient p can be assumed to be sufficiently small. In consequence, Kj is 
Lipschitz continuous. We can derive further regularity on (u,v) by (I4.6aj) - (14.6bp . and one 
concludes by a bootstrapping argument. 

Finally, the condition Kj(0, U-) = 1 follows from ( 14. 6 eft . Note here that the second 
(integral) term is of order 0(s) \Xj(u(s, s; U-)) — (t(s)\ since u', v' and d a Xj are precisely of 
this order for t G [0, s\. □ 

Remark 4.1. 1. Dal Maso, LeFloch, and Murat's definition allows for more general jump 
relations that need not be related to a specific regularized model associated with (14.11) . In 
view of the discussion above, especially Proposition \4-l\ it seems natural to fix a vector 
valued map k,j = Kj(s; uJ) and a function ) satisfying (7tj(0;tt_) = 1) and to 

prescribe the Hugoniot curve via the following differential system 

d s Hj(s;u-) = r i (W J -(s;«_)) + « i (s;«_)(A i (W j (s;w_)) - Aj(s;u_)). 
s d s Xj(s; uJ) = Kj(s; uJ) (Xj(TCj(s; — Xj(s; «_)) , 

with the initial conditions 

W i (0;u_)=u_, Aj(0;u_) = Aj(u_). (4.11) 

A^ofe the system is singular at s = 0. is however not difficult to check that (14.10j) - (14.11j) 
single out a unique (generalized) Hugoniot curve. 

2. To close this section, we point out that, more generally, it would be interesting to 
identify, within the DLM framework, suitable conditions on the family of paths ensuring 
that the Riemann problem admits a solution that enjoys the properties exhibited in the 
present paper. Recall that, for genuinely nonlinear systems, it is known that indeed the 
Riemann problem has a unique solution [27[ Wo)/ and that the Glimm scheme converges 

m 



5 Existence theorem and approximation scheme 
5.1 Existence result 

In the present section we consider the Cauchy problem for a conservative or nonconserva- 
tive, strictly hyperbolic system 

d t u + A(u) d x u = 0, u = u(t,x) G R N , (5.1) 

with 

u(0,x) = u (x), (5.2) 

where the initial data u$ : R — > W N has sufficiently small, total variation TV(uq). All wave 
speeds under consideration will remain within disjoint intervals (A™ n , A™ ax ). We establish 
here the convergence of a front tracking scheme for the approximation of (15. lft - (15. 21) . Our 
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original contribution in this section is the introduction of a wave splitting strategy and the 
derivation of new interaction estimates on wave speeds. In turn, this provides us with a 
new existence theorem. 

Following Dafermos, DiPerna, Bressan, and Risebro we construct approximate 
solutions to the Cauchy problem (15. lh - ( 15721) . which are piecewise constant and consist of 
finitely many propagating fronts. In addition to the j-fronts (1 < j < N) associated 
with one of the wave family of the system we will introduce artificial fronts of small 
total strength. More precisely, given e > we are going to construct a piecewise constant 
approximate solution u e = u € (t,x) satisfying the following properties: 

• The function u £ (t, ■) admits a finite number of discontinuities for each time t, and 
the fronts meet at finitely many interaction points, 

• The propagating discontinuities in u e are of two types: 

— j-fronts u + ) associated with a family j G {1, . . . , N} and such that u + = 
ipj(s; u_) for some s, and 

— artificial fronts («__,«+) propagating with a fixed speed Xn+i larger than A™ ax . 
We sometimes refer to such a front as an (N + l)-front. No condition is imposed 
on the jump, and the strength of such a front is defined as 

€n+i(u~,u + ) := \u+ — U-\. 

• The inner speed variation of or each j-front u + ) is uniformly small: 

$j(u-,u+) = 0(e). 



The speed of any j-front («_,«+) is uniformly close to its correct speed, that is, 
calling A the speed of the front, 

X = X j {0;s;u-) + O(e). 



The total strength of artificial fronts is uniformly small: 

\u + -u-\ = 0(e). 



artificial fronts 



We will refer to a sequence of functions u e = u e (t,x) satisfying the above properties as 
an e-approximate front tracking solution. 

Theorem 5.1. Consider a general, nonlinear, strictly hyperbolic system in conservative or 
nonconservative form (15. ip . Then there exists a constant c > such that for every initial 
data of bounded variation uq with TV(uq) < c there exists a sequence of e-approximate front 
tracking solutions which, as e — > 0, converges to the entropy solution of (I5.ip - (l5.2p . 
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5.2 Exact and approximate Riemann solvers 

We summarize the properties of the Riemann solver that were established earlier via two 
different techniques (Iguchi-LeFloch's explicit construction, Bianchini-Bressan's vanishing 
viscosity approach). The Riemann solver is an application which associates, to any two 
states Ui,u r G B§ 1 (for some small 5\ > 0), the self-similar solution u = u(£) (£ = x/t) of 
the corresponding Riemann problem. We distinguish between the three sets of conditions 
which we now describe. 

Wave curves. 

• The Riemann solution is determined from N wave curves ipj = ipj(s; u) (1 < j < N), 
which are globally Lipschitz continuous and locally different iable at the origin s — 0, 
with 

ipj(0; u) = u, d s ipj(0; u) = rj(u). 

• It consists of N + 1 states Uj separated by N wave fans (uj, %+i), determined by 

u r = ^ n {sn] ^n-i{sn-i, ■ ■ , i>i(si, Ui)- ■■ )), 
u :=Ui, Uj+i :=ipj(sj;Uj), u N :=u r . 

The parameter values Sj =: 6j(uj,Uj + i) = 6j(ui,u r ) are refered to as the wave 
strengths, and s = (s\, ■ ■ ■ , Sjv) is refered to as the strength vector. We also use 
the notation 

\P(s; ui) := i/) N (sn) ^jv-i(sjv-i, • ' ' , Ui) • • • ))■ 

• Each j-wave packet (uj,Uj + i) = (uj,i/jj(sj;Uj)) consists of (finitely or countably) 
many shock and rarefaction waves, defined as follows. Suppose for definiteness that 
Sj > 0. The range of the Riemann solution is a subset of the wave curve s G [0, sj] ^ 
ipj{s; Uj), determined by a continuous, non-decreasing wave speed 

Its generalized inverse £ h- > (Xj) 1 (£, Sj]Uj) is a (possibly discontinuous) function 
with bounded variation (which could be normalized to be right- or left-continuous). 
The Riemann solution is given by 

{ Uj , £<Xj(0,Sj-,Uj), 
il>j(s;uj), s = X] 1 {£,s j ;u j ), (5.3) 

U j+1 = ^j(Sj] Uj), £ > Xj(Sj, Sj] Uj). 

• Furthermore, the following superposition/decomposition property holds: given 

any speed £ G (Aj(0, Sj] Uj), Xj(sj, Sj] Uj)) and introducing s' := Xj (£,Sj]Uj), the 
wave packet (uj,ipj(sj]Uj)) can be obtained by simply patching together the wave 
packet (uj,ipj(s';Uj)) and the wave packet (ipj(s';Uj),ipj(sj;Uj)). 
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Interaction estimates on wave strengths. 

To any two wave packets (u-,u + ) with u + = tpj(s;U-) and (v-,v+) with v + = ^fc(r;v_), 
we can associate the potential of interaction and the amount of interaction 

Q(u-, u + ; V-, v+) > 0, I(u-, u + ; V-, v+) > 0. 

In turn, the potential of interaction between two Riemann problems (ui, u r ) = (u , • ■ ■ , mat) 
and (vi, v r ) = (v , ■■■ , v N ) is given by 

Q(ui,u r ;vi,v r ) := ^ Q(ttj-i, Uj\ u 3 --i, Vj), 

i>3 

I(ui,u r ;vi,v r ) := /(ttj-i, Uj) Vj-i, Vj), 

and the interaction potential for a piecewise constant function is obtained by summing 
Q(ui, u r ; vi, v r ) over all the discontinuities. 

The functions Q and I should satisfy the following properties: 

1. 

Q(ui,u i+ i;vj,v j+1 ) = \ei(ui,u i+ i)\ \ej(uj, u j+1 )\, i ^ j, 
Qiuij Ui+^v^Vj+x) < \ei(ui,u i+1 )\ \ej(uj, u j+1 )\, 
Q(u ,ijj i (s,u );v h v r ) = Q(u' ,tpi(s',u' );v h v r ) + O(l)(\u - u' \ + \s- s'|), 

6j(ui,U r ) = 6j(u h U m ) +6j(u m ,U r ) +0(l)li mr , 

where k mr := I{u h u m \ u m , u r ). 

2. For u a piecewise constant function u with small total variation, we have for some 
c> 0: 

Q(u) < Q(u) - Climr, 

where v! is obtained by replacing the Riemann problems (ui,u m ), (u m ,u r ) with the 
Riemann problem (ui,u r ). 

3. Q has the following behavior with respect to a decomposition of waves: if the Rie- 
mann problem (u ,u 2 ) can be decomposed into (u , u\) and (ui,u 2 ) according to the 
property above, then the interaction functional remains unchanged when splitting 
(u , u 2 ) as follows: 

Q(u ,u 2 ;vi,v r ) = Q(u , ui,vi,v r ) + Q{u x , u 2 ;v h v r ), 
Q(u ,u 2 ;u ,u 2 ) = Q(uo,Ui, u ,Ux) + Q{u x ,u 2 \ u^u 2 ) 

+Q(u , Ui, Mi, U 2 ) + Q(Ux, U 2 ] M , Ml). 
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Since our construction requires the introduction of artificial fronts, it is necessary to specify 
how the interaction amount Q(u_, u + ; i>_, v + ) is extended to artificial fronts. When at least 
one of the two fronts is artificial, we define 

\u+ — U-\ €j(v-, v + ), artificial front (iz_, u + ) and j-front (v-, v + ), 
0, (v-, v+) artificial. 

This definition is natural since artificial fronts meet all j-waves on their right-hand side, 
but do not meet other artificial fronts nor waves on their left-hand side. 

Interaction estimates on wave speeds. 

To each wave packet (uj,Uj + x) = (uj,ipj(sj;uj)) we associate its minimum and maximum 
speeds 

Xj m (uj, uj+i) := Aj(0, Sf, uj), Xj^iuj, u j+1 ) := Xj(s j: sf Uj ), 
respectively, as well as its inner speed variation 

^( Uj ,u j+1 ) := \f B3L ( Uj , Uj+1 ) - Xf^( Uj ,u j+1 ). 

More generally, when (uj,Uj + i) is the j'-wave packet in the Riemann solution (ui,u r ), we 
use the notation 

t— min , . t— min . . -r-max / \ -r-max , . 

Xj [ui,u r ) := Xj {Uj,Uj+i), Xj {ui,u r ) := Xj {Uj,Uj +1 ), 
■&j{u h u r ) := •& j (u j ,u j+1 ). 

• When u r = ipj(s; u{) and u' r = ipj(s'] ufi, 

\${u h u r )-${u' h u' r )\ < | s '- s | + | u ;-„,|, (5.4) 
so that, in particular, r &{u,ipi(s,u)) < \s\. 

• When (u m , u r ) is a j-wave packet with u T := i()j(s;u m ), interacting with a Riemann 
solution (m, Um) then, setting u' m := ipj(s; Ui), 

•&j{u h u' m ) < •&j{u m ,u r ) + 0(1) \u x - u m \ \s\. (5.5) 

The same inequality holds if the j-wave is located on the right. 

• When Ui, u m , u r is an interaction between two wave packets (of the same wave family 
j) only: 

- If ej(u h u m ) ej(u m ,u r ) > 0, then 

•&j(u h u r ) < max(i?j(uj,u m ),i?j(u m ,u r )) 

+ (Xj m (u m , Ur) - Xj m (u h Um)) + 0(1) I(u h u m ; u m , u r ). 
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- If ej(ui,u m )ej(u m ,u r ) < 0, then 

i3j(u h u r ) <max(i? j (K i ,M m ),i? i (M m ,u r )) 

+ 0(1) min(|w m - u t \, \u r - u m \). 

Remark 5.1. The Riemann solvers described in Sections^ and\Q satisfy the above condi- 
tions with, in the first case, I = Q and, in the second case, distinct values Q and X. 

Approximate Riemann solvers. 

We now introduce several approximate solvers which will be needed at each wave interac- 
tion points. The approximate Riemann solutions are piecewise constant functions which 
we construct in two steps. First, we introduce the intermediate states, that is, the u- 
components of the waves, and, second, we specify the speeds of each waves. Note that 
the solutions contain j-waves propagating with a speed close to one of the characteristic 
speeds of the system, as well as artificial waves propagating at the (large, constant) speed 
Ajv+i- First of all, we define the intermediate states, as follows. 

• Accurate solver. The intermediate states here are determined straightforwardly 
from the ones in the exact Riemann solution (ui, u r ). 

• Approximate zj-solver for different families j > i. Given a j-front and an i-front, 

u m = "0j( s i; u i), u r = ipi(s 2 ;u m ), 

the approximate solver consists of the z-wave (ui,Ui), the j-wave (ui,u 2 ), and the 
artificial front (u2,u r ) determined by 

u :=u h ui := ipi(s 2 ;ui) : u 2 := ^j{s\] u{), u 3 = u r . 

• Approximate ii-solver for a single family %. Given three states ui, u m , u r connected 
by 2-fronts, 

the approximate Riemann solution consists of the z-wave (ui,u r ) and the artificial 
front (u r , u r ) determined by 

u r := ipi(si + s 2 ; ui). 

• Artificial wave solver. When the left-hand front is an artificial front and the 
right-hand one is an z'-front satisfying 

U r = 1pi(s] U m ), 

we introduce u m := ipi(s;ui) and the approximate Riemann solution consists of the 
i-wave (ui,u m ) and the artificial front (u m ,u r ). 



52 



The second part of our construction consists of prescribing the speeds of propagation of 
the waves. Fronts having a large inner speed variation will need to be split in two or more 
fronts. Fix a threshold e > 0. First of all, note that artificial fronts are never split and 
always travel at the speed Xn+i- To handle j-wave packets («_, u + ) = (u_, if}j(sj] «-)), we 
distinguish between two strategies: 

• No-splitting strategy. In this case, we propagate (u-, u + ) as a single front traveling 
at the smallest speed of the associated wave fan. 

• Splitting strategy. In this case, we replace the front by several smaller fronts 
defined as follows. Using the wave speed function Xj(s') := Xj(s / , Sj] uJ) we introduce 

p ._ - A,-(0) + 1 

■= Aj(0) + p-iMsj) ~ A;(0)], V = 0, . ■ ■ , P - 1, 

where [a\ denotes the greatest integer less than or equal to a > 0. Then, along the 
wave curve from it_, we pick up the states associated with the speeds /x p , that is 

w p := ^(sj,;«_) 

min{s' G [0, Sj], Xj(s') = /x p }, s y > 0, 
maxjs' G [sj,0], Ay(s') = < 0. 
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Then the z-wave packet is approximated by P fronts with small strength, as follows: 

fU_, x/t<fL , 
w P , fip-i < x/t < fi p , (p= 1,...,P- 1), 




where, for simplicity in the notation, the Riemann solution has been centered at the 
point (t, x) = (0,0). Note that each front corresponds to an exact solution of the 
system of conservation laws (in the sense that the right-hand state lies on the wave 
curve issuing from the left-hand state) although its propagation speed (in general) 
only approaches the true speed. 

5.3 Front tracking approximations 

We are now in a position to describe our algorithm, and we fix e, 5 > satisfying 

5 « e 2 . 

Let Uq be a piecewise constant approximation of u such that 

H - "olU 1 < e, TV ( u o) < TV(u ), (5.8) 
u s q contains 1/e discontinuity points, at most. (5.9) 
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At each discontinuity point of u e we use the accurate solver described earlier together 
with the splitting strategy and we define the approximate solution u e = u e (t,x) locally in 
time. The solution is extended until two fronts meet. As usual, the speeds of the fronts may 
need to be modified (cf. Remark 15.21 below) in order to ensure that only two fronts meet 
at every interaction. To extend the approximate solutions, we distinguish between several 
types of interactions and we use either the accurate or the approximate solver together 
with the splitting or the no-splitting strategy. 

As a rule, only fronts with inner speed variation less than or equal to e will be generated 
by the scheme. However, after one or several interactions, the inner speed variation may 
have increased so much that it is greater than 2e and need to be split in two fronts. 

• Large interactions. We call large interaction an interaction that involves an i-front 
(ui,u m ) and a j-front (u m ,u r ) such that 

I(ui,u m ;u m ,u r ) > 5. 

The solution is defined beyond the interaction time by using the accurate Riemann 
solver, together with the following rule for the speeds of the fronts : 

— For the outgoing i- or j-waves that have an inner speed variation < 2e, we 
use the no-splitting strategy. 

— For all other outgoing waves we use the splitting strategy. 

• Small interactions. We call small interaction an interaction that involve an i-front 
(ui,u m ) and a j-front {u m ,u r ) such that 

I(u h u m ;u m ,u r ) < 6. 

The solution is defined beyond the interaction time by using the approximate ij- or 
ii-solvers. The speeds of the fronts are determined exactly as in the case of a large 
interaction (recall that artificial fronts are never split). 

• Artificial wave interactions. At an interaction between a left-hand artificial front 
and a right-hand j-wave, we use the artificial interaction solver, with again the same 
splitting/no-splitting strategy on the j'-th outgoing wave. 

Remark 5.2. The speeds of certain fronts may need to be slightly modified in order to avoid 
both interaction points involving more than two fronts and interaction times involving more 
than one interaction. Starting from t = or a time t at which an interaction has occured, 
we let the fronts evolve at the speed determined earlier. We consider the first time r > t 
where two fronts meet. If there are only two fronts involved, there is nothing to do. If 
there are three (or more) fronts interacting at the time r (possibly at different locations), 
we consider the most-left fronts a and (3 (with a on the left of (3). Now, if a is a j- 
front (1 < j < N), we increase the speed of a by an amount less than e, so that the 
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new interaction time is beyond the time t. If a is an artificial front, then (3 must be a j- 
front, and we diminish the speed of f3 by an amount less than e so that again that the new 
interaction time is beyond time t. By this procedure, in both cases we have only decreased 
the interaction time. Hence, as r was the first interaction time greater than t, the new 
interaction time t' G (t, r) corresponds to a single interaction of two fronts. An important 
consequence of our construction is that two fronts that meet in the process would have met 
regardless of the speed modification. 

Based on the properties of the Riemann solver listed at the beginning of the present 
section, it is a standard matter to derive uniform estimates for the following Glimm-type 
functionals: 

V(t) :=^|e a |, Qit) :=^2Qap, 

a a,P 

where the summations are over all fronts a, (3 in the approximate solution at the time t. 
Here, e a denotes the strength of the front a and Q a p is the interaction amount between 
two fronts a, (3. 

Proposition 5.1. There exist constants u, c,Cq > so that for every initial data satisfying 
TV(uo) < v and for every sequence Uq satisfying (I5.8p - fl5.9p . (as long as the number of 
waves and the number of interaction remain finite) the quantity (V + CoQ)(t) is non- 
increasing in t and, more precisely, at every interaction involving some fronts (ui,u m ), 
(u m , it r ) 

[V + C Q](t) 

:= (V + C Q)(t+)-{V + C Q)(t-) 

(5.10) 

■c rmn(\u m — Ui\, \u r — u m \), non-monotone interactions, 
■cl(ui,u m ;u mj u r ), otherwise. 



< 



5.4 Total number of fronts and interaction points 



In view of Proposition 15. 1[ to establish that the algorithm is well-behaved we need only 
show that the numbers of fronts and interaction points remain finite. Define 



N 



Af(t) := ^A4(t), Afj{t) := #{ j-fronts at time t}, 

k=i 

N [•& -e) 

e £ (t) : = 5> £)fc (t), e ed (t) := £ K a £ >\ 

k=l j'-fronts 

and introduce the functional 

F{t) := d (V(t) + Co Q(t)) + 3 £ + JV(t), 
where C\ is some positive constant and Cq is the constant already introduced in (15.101) . 
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Lemma 5.1. For 5 sufficiently small and C\ large enough (defending on e and 5) the 
functional T is non-increasing. 

Proof. The function T is constant away from interaction times; we distinguish between 
every type of interactions. Each interaction involves a left-hand, z-front a = (ui, u m ) with 
strength s a (1 < i < N + 1) and a right-hand, j-front j3 = (u m ,u r ) with strength s@ 
(1 < j < N). We denote by 9 a ,X a and 6p,\^ the inner speed variation and wave speed 
function associated with the incoming waves. Call 77^ the outgoing k- waves (1 < k < N+l), 
with inner speed variation d Vk and wave speed function . Finally, call X the following 
modified amount of wave interaction: 

• X :— I (the interaction amount at the point), unless the interaction is a non-monotone 
interaction of waves of the same family, 

• X = the strength of the smallest incoming wave, if the interaction is a non-monotone 
interaction of waves of the same family. 

We will summarize the evolution of the terms arising in T in a table (see Table [1]), 
according to the nature of the interaction, the family of the wave considered, and whether 
the splitting strategy is applied or not. In this table, "z 7^ j" (respectively u i = j") 
refers to interactions between waves of different (resp. same) family; "acc" (respectively 
"app.") stands for "accurate" (resp. "approximate"). We regard artificial interactions as 
approximate interactions between waves of different families. Since artificial fronts are not 
considered in the last two terms of J 7 , they are not included in the table. 

We fill up the cells of the table by relying on our interaction estimates on the strengths 
and the inner speed variation, and the definition of the number of fronts when the splitting 
strategy is applied. The following remarks are in order. 

1. In the case of a monotone interaction between two fronts (ui,u m ) and (u m ,u r ), 
the term (A™ m (w m , u r ) — \™ m (ui, u m )) + in (15.61) vanishes. This is due to the fact that 
fronts evolve at the minimum speed of the wave packet - except in the case described in 
Remark 15 .21 In the latter case, however, two fronts that meet would have met regardless 
of the slight change in speed. 

2. When the i-th or j-th outgoing wave is split, we have (for instance for the z-wave 

a): 

• In the case of a monotone interaction of waves of the same family, 

IW)] = L^J - 1 

• In the case of a non-monotone interaction of waves of the same family, 

[W)] = L^J - 1 

< -1 + — + 0{-)I(t), 
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Otherwise, 



Moreover, we have 



Wt)] = L^J<7 + o(iw 



„ — ( maxima, 'da) — e)/e, non- monotone case, 



The fact that the wave was split implies that •Q rii > 2e. From the estimate 



>h — 



max(i9 a , ftp) + 0(l)X(t), non-monotone case, 
tf a + 0(l)J(t), otherwise, 



we deduce, in both cases, 

[e £ii (t)] <-i + o(I)j(t). 

Using again (15.111) . this yields 

3Te < -1 + 0(-)I- | max (^"'^)/ £ ' non-monotone case, 
~~ £ \'d a /£, otherwise. 

Now, the desired conclusion follows from the inequalities summarized in Table [U At 
each interaction, we have [V + CqQ~\ (t) < — (Co/2)X(t), as follows from Proposition 15.11 
Hence, it is sufficient to take C\ = 0(l/e), except for the N — 2 (or N — 1) contributions 
of 1 in Af(t) due to new fronts of families k ^ i,j in the accurate solver. But, such fronts 
appear only with the accurate solver; hence, when X(t) > 5, it suffices to take C\ = 0(1/ S) 
to cover all cases. □ 

Number of fronts and interaction points. 

• Finite number of j -fronts (1 < j < N). It follows immediately from Lemma ISTTl that 
the total number of j-fronts is bounded for all times. Moreover, observe that at all 
interactions generating new j-fronts an amount > min(l, Co Ci5/2) is used out of 
J-{t)] this implies 

M(t) < i. 

• Number of accurate interaction points. From Proposition 15.11 it follows also that the 
number of accurate interactions is finite. 

• Number of interaction points using the splitting strategy. At such points, as seen in 
the table, T decreases of an amount of 1 at least; hence, this case arises finitely many 
times. 
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Interaction 


Wave /c 


Strategy 




3[e e , fc ] 


% ^ j, acc. 


f or j 


CT~\ 1 1 'j - 'j - 1 (T 

dJJIIL Lilllg, 




< i 

— e 

+ 0(l/e)I 


no splitting 





< 0(l/e)T 




splitting 


1 + 0(l/e)Z 





i — j, acc. 


i 


GTll l"t""t"lTl CT 
SJJ11L Llllg, 


^ max(() a 

- 1 + 0(l/e)l 


— e 

- 1 + <D(l/e)I 


no splitting 


-1 


< 0(l/e)X 


k j£ i 


splitting 


1 + 0{l/e)I 





i 7^ j, app. 


i or j 


QTH 1 1~ t~ 1 T 1 O' 
BLJ11L Llllti 




< - 1 

— e 

+ C?(l/e) J 


no splitting 





< 0(l/e)I 




splitting 








» = j\ app. 


i 


mil ittiTi p 1 


^ max(() 
— e 

- 1 + 0(l/e)l 


^ max(i9 a ,#^) 
— e 

- 1 + 0(l/e)J 


no splitting 


-1 


< 0(l/e)T 


k ^ i 


splitting 









Table 1: Estimates across each type of interaction 



• Number of approximate i-interaction points using the no-splitting strategy. This num- 
ber is also finite since a z-front is lost at these points. 

• Number of approximate ij interaction points using the no-splitting strategy. When 
tracing forward a given front by taking any front of the same family as its successor, 
we see that two fronts of differents families that have met will not meet again, since 
the speeds of different families are separated. Since interaction points from which at 
least three j-fronts leave, are in finite number, this yields that approximate ij points 
and artificial interaction points are finite in number. 

• Number of artificial fronts. There are only a finite number of approximate interac- 
tions and, therefore, the total number of artificial fronts is finite. 

• Number of artificial interaction points using the no-splitting strategy. The same con- 
sideration as for the approximate ij interaction points using the no-splitting strategy 
applies. 
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5.5 Number of fronts of early generations. 



We already have a bound on the total number of fronts, but it will be crucial in the 
next paragraph to have a sharp bound on their number, if we restrict to fronts of early 
generations, as we describe now. 

Front generations. 

We define the generation g a of a front a as follows. 

- All fronts outgoing from the initial line have generation 1. 

- At an zj-interaction (an z-front a meeting a j-front j3), the generation of outgoing 
i-fronts (resp. j, k ^ i,j) is fixed as g a (resp. gp, max(g Q , gp) + 1), whether or not 
these fronts are split. 

- At an z-interaction (of two z-fronts a and j3), the generation of the outgoing i- fronts 
(resp. k 7^ i) is fixed as min.(g a ,gp) (resp. max(g a , gp) + 1), whether or not these 
fronts are split. 

- At an artificial interaction point the generation of each family is preserved. 
Main estimate. 



Two fronts are said to be approaching if they are of the same family, or else if the front of 
the largest family is on the left. Hence, V k represents the number of potential interactions 
between fronts, with largest generation number equal to k. Note that artificial fronts 
are taken into account in the above quantities (except for what concerns the inner speed 
variation). Since there is no front of generation 0, all these quantities vanish for k — 0. 
Estimates on the number of fronts of generation < k are now derived. 

Lemma 5.2. There exist a constant Ci = (^(e) > ; an increasing sequence A k = A k (e) > 
0, and two decreasing sequences Bk = Bk(e) > and Dk = D^{e) > 0, such that for all 
k > the functional 



We define 



AT(r) := #{ fronts at time r, of generation fcj, 
^P k '■= #{( a i(3) fronts at time r/ max.(g a , gp) = k, a and (3 approaching}, 




a front of 
generation <k 



k fc-1 



T k := C 2 A k (V + C Q) + 3A k Bf k + ^ BjT + ^ D-P* 




i=l i=l 
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is non-increasing in time. This is true regardless of the value of 5 (small enough with 
respect to e). Moreover, for each k, A k (e) is bounded above by a polynomial expression in 
1/e, while B k (e) and D k (e) are bounded below by polynomial expressions in e. 

Corollary 5.1. For each k, the total number of fronts of generation less than or equal to 
k is 

k 

^AT<J fc (e), (5.13) 

i=l 

where Jk(s) is polynomial in 1/e. 

Proof of Lemma 15.21 and Corollary 15. 1L 

We start with two observations. First, from the proof in Section that for some C 2 = 
C 2 (e) = 0(l/e), the functional 

g k :=c 2 (v + c Q) + 3ef\ 

is non-increasing in time. Hence, if the result in the lemma is established for a given A^, 
then it remains true with a larger A^. Moreover, as seen in Lemma \5. 11 Q k decreases by 1, 
at least, when an incoming wave is split. 

Next, when an interaction creates new fronts, the total number of new fronts is at most 
0(1/ e), as follows from (15.41) and the fact that the total strength of the waves is bounded. 

Now, the proof is done by induction on k. The case where k — 1 is essentially given by 
Lemma [5.1 1 In fact, when the fronts of generation g > 2 have not to be considered, the 
estimate on C\ can be replaced by C 2 = 0(1/ e). (Note that there is no artificial front of 
generation 1.) We can choose for instance D\ — 1 and B = 1. 

From now on, let us only consider the passage from k to k + 1, which is the heart of 
the proof. Suppose that for a certain k > 1, the expression in (I5.12p is non-increasing. 
Whenever the property stated in Lemma \5. 21 is established at the rank k, Corollary 15. II at 
the rank k follows immediately; hence, we can use ( 15. 13ft at the rank k. 

The goal is to determine A k +i (large enough) and B k +i, D k (small enough) so that 
the desired property is valid. We distinguish between several cases when one of the "new" 
terms (that appear in J-~k+i but did not appear in jF fc ), that is, J\f k+1 and V k , can grow. 
This will provide us with some conditions on A k+ x, B k+ i and D k , that we will be able to 
fulfill by choosing these constants sufficiently large or small. 

Consider the interaction of two fronts with corresponding amount of interaction /. We 
introduce the same notation X as previously (including the case of artificial interaction), 
that is, X := Q (the potential at the interaction) unless the interaction is a non- monotone 
interaction of waves of the same family, in which case X is the strength of the smallest 
wave. There are two cases where V k can grow, and two cases where J\f k+1 can grow. 

1. Cases that increase V k . 
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The value V k can grow in two ways: (1) an existing front of generation < k is split after 
an interaction with another wave, or (2) a new front of generation < k is created after an 
accurate interaction where the maximal generation of the incoming waves is exactly k — 1. 
Note that we can measure the increase of V k by 

[P k ] < (number of new fronts ) x J k {e). 

• First case: splitting of a front of generation k. In that case, the number of new fronts 
is of order 0(l/e); we can compensate the increase of V k by taking A k > J k (e)/e, 
since the decrease of Q k in that case is at least 1 . 

• Second case: creation of a new front generation k. We suppose that no incoming 
wave is split (if not, enhancing A k allows us to absorb this increase too.) Then, V k ~ x 
decreases by 1, at least, and it suffices to consider D k < eD k _i/ J k (s). 



2. Cases that increase AT 

The value Jx can grow in two ways: (1) an existing front of generation k + 1 is split 
after an interaction with another wave, or (2) a new front of generation k + 1 is created after 
an accurate interaction where the maximal generation of the incoming waves is exactly k. 

• First case: splitting of a front of generation k + 1. In that case, the increase of the 
term N due to these new fronts is compensated by the decrease of A k+1 Q k+ i (at 
least A k+ i), provided A k+1 > B k+1 as seen in Tabled! 

• Second case: creation of a new front generation k + 1. Again we consider only the 
case where there is no splitting of incoming waves. The increase of the term Jx 
due to these new fronts is of order A^ — 2 + 0{l/e)X {ox N - 1 + 0(l/e)T). In that 
case, J\f k (if incoming waves are of same family) or V k (otherwise) decrease by 1. 
Hence, taking B k+ i small with respect to B k and D k allows to get the decrease. 

This concludes the proof of Lemma 15.21 and Corollary 15.11 □ 
5.6 Conclusion 

The fact that the scheme converges to the entropy solution in the limit is a consequence 
of the following two properties: 

• The j-fronts (1 < j < N) travel approximately at the correct speed given by the 
Rankine-Hugoniot relation. 

• The total strength of artificial fronts remains uniformly small. 
Both properties are now discussed. 
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Accuracy of the speed of j-fronts. 

From our construction it follows that, for any j-front a, 

d a < 2e. 

Indeed, we have: 

• When the front under consideration is generated by the accurate solver (or by the 
initial solver), then d a < e. 

• When the front is generated by an approximate solver, then either $ a < e (if the 
solver happens to split the wave of the family of a) or d a < 2e (otherwise). 

Furthermore, the front travels at a speed which is the lowest speed in the wave packet up 
to a e error at most. So, calling A the speed of the front and a a (s) any speed in the wave 
packet, we also obtain 

\X-a a (s)\ < 3e. 
Total strength of artificial fronts. 

We follow here the argument known in the genuinely non-linear/linearly degenerate case 
(see, for instance, [10], Section 7.3.6). 

Strength of an artificial front. From our previous discussion it follows that, for any artificial 
front a, 

\e N+1 (a)\=0(l)5. 

Indeed, with the previous convention and thanks to the interaction estimates, the strength 
of any new artificial front is of order 5. One can follow it during successive interactions. 
Calling V a the total strength of fronts approaching a (that is, all j-fronts on its right), one 
gets (PUj, p. 139) 

\S a (t)\<O(l)5exp(C'(V a + C Q)). 

Total strength of artificial fronts. For k > 1, we call (resp. V k art ) the total strength of 
fronts (resp. of artificial fronts) of generation > k, and Qk defined as previously, where the 
sum is over all couple of fronts a and (3 for which max(g Q , g@) > k. 

The estimates on the strengths of waves according to their generation remain valid; 
hence, provided the total variation is small enough, we can deduce that, for some 7 < 1 
and for all times t, 

Q k (t) < c 3l k , v k (t) < c a k - (5.14) 

Now, given any integer M we have 

v art (t) = v k rt + Yl y k rt - 

k<M k>M 
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The first term is estimated by 



J2v k art <0(l)5J M (e), 



k<M 



while, for the second one, 



Af+l 



£^<C4f- 

k>M ' 

Hence, given e > 0, we can choose M so large that the second term above is less than e/2, 
and next we choose 5 small enough so that the first term is less than e/2 too; hence 



V art (t) < e. (5.15) 



Convergence. 



Using the L%°(BV X ) bound and the uniform bound on the wave speeds, it is a standard 
matter to derive a Cip t (L x ) bound. Relying on Helly's theorem, these estimates allow us 
to extract a converging subsequence, say 



u^ue Ll c (R N ) 



We now check that u is an entropy solution. We will consider here the case of a conservative 
system endowed with a convex entropy pair and prove that for any non-negative test- 
function if : x K — > R and for every smooth, convex entropy / entropy-flux pair 
(77, q):R N ->Rx R N . 

POO p 

liminf / / (r](u € ) tp t + q(u e ) <p x ) dtdx > 0. (5.16) 

™^+°° Jt=0 JxGR 

Calling /„ the above integral, for supp(</?) C [0, T] x IS we have 

In = I i^Mt) - [q(u e )} a (t))^(t, x a ) dt, 

where x a is the trajectory of the front a and [h] a is the jump of a function h on the front 
a. Thanks to the regularity of rj and (15.151) . the sum over artificial fronts is 0(e). The 
statement (15.161) follows immediately from: 



Lemma 5.3. For any j-front a and all times T\ and T2 one has 

/•T2 



> 



M*) [V(u £ )]a{t) ~ [q(u e )] a (t)Mt,x a )dt 

(5.17) 
-e\s a \\T 2 - Ti| II^Hci- 
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Proof. Consider a front a connecting w_ to u + = tf}i(s a ,U-), and denote by Aj(-,s Q ,M_) 
the corresponding wave speed function. Recall that the corresponding solution is given in 
(15.31) . Denote by lj the corresponding Riemann solution centered at the point (Ti,x a (Ti)). 
The solution satisfies 



Ti 



(r](uj)(f t + q(uj)ip x ) dtdx 
We also introduce, for t G [Ti, T 2 ], 

co a (t,x) = 



rj(uj)ip(-, x) dx 



Tz 



> 0. 



(5.18) 



m_, cc/t < A c 
M + , x/t > A c 



where \ a is the speed of the front a (that is, X a = A,,(0, s a , w„) up to the additional change 
of speeds in Remark |5~21) . Clearly the left-hand side of (I5.17P is equal to the left-hand side 
of (15.181) when we replace Co by u a . Hence, it is sufficient to prove that each term in the 
the left-hand side of (15.181) yields an error of order £|s a ||T2 — Ti | H^Hc 1 when one replaces 
lj by uj a . 

For the first integral, the difference is supported in a triangle whose area is of order 
0{l)e{T 2 — Ti), and the difference of the integrand is of order 0(1) |s a |. For the second 
integrals, the difference is supported in an interval whose length is of order 0(1) e(T2 — Ti), 
and the difference of the integrand is again of order 0(1) \s a \, which concludes the proof 
of Lemma 15.31 □ 



6 Time- regularity of graph solutions 

We now use the front tracking scheme to study the regularity of graph solutions introduced 
in [30]. To simplify the presentation it is convenient to assume that the flux f(u) is defined 
for all u G M. N . To begin with, we need a few definitions from [30]. (Note that Lipschitz 
continuous representatives of all geometric maps under consideration are used throughout 
the present section.) 

6.1 Geometric version of the front tracking scheme 

A parametrized graph is a map (X, U) : M — > R x such that X and U are Lipschitz 
continuous, and 

d s X > 0, lim X(s) = ±oo. 

s— >±oo 

A maximal interval [s_,s+] where X remains constant will be refered to as a vertical 
segment of the graph. We say that the parametrized graph (X, U) contains a single 
shock if the associated BV function U o X^ 1 is a step function with a single discontinuity 
point. A (Lipschitz continuous), time-dependent parametrized graph is a map (X,U) : 
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I + xl^lxR w such that X, U G L°°(R+, £ip{R)), (A, U){t) is a parametrized graph in 
the sense above, and for every continuously differentiable function 9 with compact support 
and uniformly with respect to t > 0, 

/ Ud t 9(t,X)d s Xds < ||0(V)||l~(k). (6.1) 

A DLM family of paths (Dal Maso, LeFloch, and Murat [TB]) is a Lipschitz contin- 
uous map $ : [0, 1] x R N x R N such that for all ui, u r , u[, u' r G R N 

§(0;u h u r ) = u h u r ) = u r , 

\\d s ®(-;ut,u r )\\L°°{o,i) ^ \u r -ui\, (6.2) 
\\®(-;ui,u r ) - ^(■;u' h u' r )\\ L oo ( p j i) < \u r - u' r \ + \ui - u\\. 

A family of Riemann graphs is a map (A, $) : [0, 1] x R^ x R N -^Rx R^ such that : 

• For any ui,u r G R N , (A(-;ui,u r ),$(-;ui,u r )) is a parametrized graph, and $ is a 
DLM family of paths with moreover 

\\A(-;ui,u r ) - A(-;m;,mJ,)|| L oc (0i1) < \u r - u' r \ + \ui - u\\. 

• The system of conservation laws (15. ip is satisfied: for every smooth function 9 : R — ► 
R and all u h u r G R N 

/ (— A(s; itj, it r ) + Df o $(s; u;, w r )) d s $(s; u;, w r ) 6 1 o A(s; itj, u r ) = 0. (6.3) 

•/ 

• Along every vertical segment [s_, s+], the Riemann graph (A, $)(•; it-, u+) associated 
with u± := <&(s±; M r ) contains a single shock and 

$(■; i*r)|[ a _, 8+ ] = $(■; «+) o /?, (6.4) 

where (3 : [s_, s + ] — > [0, 1] is the linear map satisfying = 0, /3(s + ) = 1. 

Definition 6.1. A family of Riemann graphs (A, $) being fixed, a parameterized graph 
(A, Z7) : R+ x R -> R x R N is called a graph solution of dSU) subordinate to (A, $) 

• For ew en/ test-function 9 : R + xl->l, 

( - U d t 9(t, X) + /([/) d x 6(t, A)) <9 S A dscft = 0. (6.5) 



+ XJ1 



For almost every time t and on every vertical segment [s_,s + ] the graph solution 
coincides with a prescribed paths, more precisely 

U{t\ =$( M „, M+ )o/3, (6.6) 
where u± := U(t, s±), and (3 : [s_, s + ] — > [0, 1] is the linear map satisfying /3(s_) = 0, 
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Remark 6.1. (i) The condition (16.11) is equivalent to saying that the BV function u = 
[/or 1 , belongs to Lip(R + , L 1 (R)), that is 

\\u(t) -u(t')\\ LHR) < \t-t'\. 

(ii) The condition ( 16. 5ft is equivalent to say that u — U o X~ l is a weak solution of 
(15. ip in the sense of distributions. 

(Hi) In Definition \6.1\ the same parametrization (up to linear rescaling) is used for the 
graph solution and the prescribed paths. This definition will be sufficient for the purpose of 
the present paper. 

To implement the front tracking scheme we will need to restrict the class of Riemann 
graphs, relying here on the interaction amount Q(ui, u m , u r ) associated with three constant 
states. 

Definition 6.2. Consider the strictly hyperbolic systems of conservation laws (15.11) . A 
family of Riemann graphs (A, $) is said to satisfy the interaction estimates if for all 
ui,u m ,u r G ~R N such that u\ is connected to u m by an i-wave fan and u m is connected 
to u T by an i-wave fan 

|| (A, u r ) - (A, u m ) V (A, $) (•; u m , u r ) ||l 1 (o,i) ^ Q(ui,u m ,u r ), (6.7) 

where (A, $)(•; ui, u m ) V (A, $)(■; u m , u r ) : [0, 1] — > ~R N is the arc-length parametrization of 
the concatenation of the two maps. 

By the result in the earlier section, there exists a family of Riemann graphs satisfying 
the interaction estimates. 

The main result in this section is as follows. 

Theorem 6.1. (Geometric version of the front tracking scheme.) Let (A, $) be a family 
of Riemann graphs satisfying the interaction estimates. Let u e = u e (t,x) be a sequence 
of front tracking approximations for the Cauchy problem (l5.1l) -( !5T2l) . Then there exists 
a parametrized graph (X e ,U e ) such that X e is Lipschitz continuous in both (t,s), U e is 
continuous in the space variable s, and the following uniform estimates hold: 

\\(d t X*,d s X*)\\ Lao{R+xR) < 1, HWUc^xR) < TV(u ), 
\\U^t)-U^t')\\ L , m <\a^t)-a^t% M'>0, 1 • j 

where a e : R + — > M + is a non-increasing, uniformly bounded function measuring the amount 
of cancellation and interaction in u e . Moreover, (X e ,U e ) converges uniformly for all but 
countably many times t toward the graph solution (X, U) of the Cauchy problem (15. II) - 05. 2p . 

Furthermore, the map X is Lipschitz continuous in both variables (t, s) while U is 
Lipschitz continuous in s and satisfies 

\\U(t) - U(t')\\ LHR) < \a(t) - a(t')\, t,t' > 0, (6.9) 

where a : M. + — > R + is the pointwise limit a e . In particular, U is continuous at all but 
countably many times t. 
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Remark 6.2. It is known that the BV solution u of the problem (j5.ip - (j5.2p belongs to the 
space Lip(M + , L 1 (IR)) ; which is expressed by the condition (16. 3p . In contrast, the estimate 
( 16.91) provides a control of the component U even in the vertical segments. 

6.2 Proof of convergence 
Step 1. 

We begin by denning the parameterization. The standard approach consists of introducing 
the arc-length parametrization of the graph of u e , which is based on the total variation 
map x i — > TV r _5 00 (ti e (t)). In order to ensure that the graph be continuous in time and 
since the total variation of u e may change at times where two fronts meet, we need to 
modify the arc-length parametrization as follows. We will take into account the interaction 
and cancellation of waves, and take advantage of the decrease of the generalized Glimm 
functional. To the sequence u e , let us associate the non- negative measure fj, e consisting of 
Dirac masses concentrated at the points (t Q ,xo) where two fronts (ui,u m ), (u m ,u r ) meet. 
The point mass at (to^o) is the scalar fi e (t ,x ) given by 

//(t ,x ) := -(Q e (t +) -Q e (t -)) > 0. (6.10) 

In view of the results in previous sections the total mass of this measure is uniformly 
bounded 

»xR)<l. 

Let Aqo be a large positive constant, larger than all speeds Xj as well as the speed Aat+i 
of artificial waves introduced in the front tracking scheme. Consider the map 

a e (t,x) :=x + V e (t,x)+fi e (n t , x ), (6.11) 

where V e is the total strength of waves at time t on the left of x, and we have introduced 
the triangular region 

fit,* : = {(s,y) I <s <t, y < x- Aoo(t - s)}. 

Clearly, for every time t, the map c € (t) is strictly increasing in x. Moreover, a e is continuous 
in time except along polygonal lines which are transverse to the line t =constant. 
Let us now define X e : R + x E -> R by 

X e := (a 6 )- 1 . (6.12) 

By the properties of a e just mentioned and the fact that u € is piecewise constant, it is not 
difficult to check that X e is Lipschitz continuous in both t and s. 

Finally, for every t that is not an interaction time, following [16] we define U e (t) from 
u e (t) and X e {t) by completion based on the given family of paths $. Precisely for each t, to 
the family of paths $ and the BV function u e (t) one can associate a unique locally Lipschitz 
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continuous, parametrized graph (X e , U e )(t) called the "^-completion of u e (t), characterized 
by the conditions that 

u e (t) = U e o (X 6 )- 1 ^), (6.13) 
and, on every maximal interval [s_, s + ] on which X e (t) is constant equal to some value x, 

U\t\ = $(u e (t,x-),u e (t,x+)) o(3, (6.14) 

where (3 : [s_, s + ] — ► [0, 1] is the linear map satisfying /3(s_) = 0, /3(s+) = 1. 

The regularity in space is clear by construction, since u e has uniformly bounded total 
variation in space and we are using (a modification of) the arc-length parametrization. 
The function dtX e behaves essentially like d s X e in a neighborhood of a single wave front. 
On the other hand, the function X e remains constant at an interaction point. This yields 
easily the uniform bounds for d s X e , d t X e , d s ll e . 

To control d t U e we observe that 

\U%s)-U e (t',s)\ < \t-t'\ sup\d s U(t,s)\ 
<\t-t\supTV(u% 

provided there is no interaction point within the time interval [t, t']. Thus we only need to 
discuss the behavior of dtU e at interaction times. 

Consider a point (to,Xo) of interaction of two wave fronts, (ui,u m ) and (u m ,u r ), and 
denote by := [sj, s m ], I r = [s m , s r ] the s-intervals describing the incoming fronts, with 

Sm |6("U;, U m j | , S r S m |f(llm] ^r) | • 

In view of (16.101) and (16.111) . the interval / = [s;,s r ] is used to parameterize the outgoing 
wave fan. Let us begin with two waves of different families. Comparing the graphs before 
and after the interaction we can write 

\\U e (t +) - £/ e (to-)IU~ W = \\U e (t +) - U e (t -)\\ Loa{I) 

< min(s m - s h s r - s m ), 

since whenever the left-hand side is Lipschitz continuous in the data Ui, u m , u r and vanishes 
whenever the right-hand side vanishes. By integration in s this inequality leads us to the 
L 1 bound 

\\u*(t +) - u € (t -)\\mR) = \\u £ {t +) - u^t -)\\ LHI) 

< ( Sr - Si) \\U e (t + ) - f/ e (t -)||L-(/) 

< (s r - si) min(s m - s h s r - s rn ) 
<2(s m - si) (s r - s m ) < \u m - ui\ \u r -u m \. 

Since the interaction potential between two waves of different families is quadratic we have 
obtained 

||^o+) - U £ (t -)\\ LHR) < Q e (t -) - g e (t +). (6.15) 
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At a point where two fronts of the same family meet, we need a more precise bound 
and we rely on the property (16.71) satisfied by the family of Riemann graphs. We have 
immediately 

wu*(t +) - [/ e (t -)iiLX( K) = \mt +) - U^t -)\\ LHI) 

<Q(t -)-Q(t +), 1 ■ ] 

which completes the derivation of uniform estimates. 
Step 2. 

We now justify the passage to the limit e — > 0. On one hand, let us denote by u the 
limit of the sequence of front tracking approximations u e . On the other hand, consider the 
parametrized graphs {X e , U e ). In view of the uniform Lipschitz bound on X e (t) (see (16. 8h ) 
and after extracting a subsequence if necessary we can assume that there exists a Lipschitz 
continuous function X such that 

X e (t) — > X{t) uniformly on all compact subsets in (t, s), 

d s X%d t X e d s X,d t X weak-* in ( ' ' 

On the other hand, since U e and d s U e are uniformly bounded we can find U(t) defined for 
all rational numbers t such that 

U e (t) ->> U(t) weak-* in W^°° (6.18) 

for all rational t at least. Using the uniform bound on dtU e (see (16.81) ) it follows that U can 
be actually defined for all times and that the convergence above holds at all but countably 
many times t. Let us set w := U o X. 

For t fixed the inverse map (X e ) _1 (t) is monotone increasing and thus converges point- 
wise some some limit Y . From the identity X e o (X 6 )' 1 = id we deduce that X o Y = id, 
so that Y is the (generalized) inverse of X. In turn, passing to the limit in the relation 
U e o (X 6 )" 1 = u e yields w — U o X~ l = u. In other words the BV function associated 
with the limiting graph is exactly the function u. In particular, the graph satisfy the 
conservation law in the sense ( 16.51) . 

To conclude it remains to determine the vertical parts of (X, U) and establish (16. 6ft . We 
will check that (X, U) coincides (up to re-parametrization) with the ^-completed graph of 
u whose arc-length parameterization is denoted by (Y, V$): 

{X,U) = (Y,V*)op, (6.19) 

where (3 : K — > M is an increasing, onto, and Lipschitz continuous map. In particular, 
(I6.19P implies that the vertical parts in both graph coincide so that (16. 6p holds. 

The statement (I6.19P is a non-trivial property of the front tracking scheme which we 
will establish by relying on the well-known property of local uniform convergence : for 
all but countably many times t and for each r\ > 0, any point ifK must be 
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1. either a point of continuity of u(t) for which there exist a neighborhood N v (x) and a 
real > sufficiently small so that for all y G N v (x) and e < e v 

\u e (t,y) - u(t,x)\ + \u(t,y) - u(t,x)\ < rj, (6.20) 

2. or else a point of jump for u(t) for which there exist a neighborhood N v (x), a suf- 
ficiently small real e v > 0, and a sequence of points x e G N v (x), so that for all 
y G N v (x) and e < e,, 

|u e (t,y) - u(t,x±)| + |u(t,y) - u(t,x±)\ < f] \i y<,x\ (6.21) 

Specifically, the argument presented now shows that the local uniform convergence of 
functions implies, for the associated ^-completed graphs, the convergence in the uniform 
sense of graphs. 

Fix a time t at which u € (t) — > u(t) uniformly locally. Given t] > we can select finitely 
many points Zi, . . . , z m so that 

\u + (t,x)-u-(t,x)\<r). (6.22) 

x£Rx^zi,...,z m 

To each Zj we associate a neighborhood N v (Zj) and a point so that the property f)6.2ip 
holds at Zj with 77 replaced by rj/m: 

\u € (t,y) - u(t, Zj±)\ + \u(t,y) - u(t, Zj ±)\ <r]/m 

if 3/^,2/ eN v ( Zj ). [ ' ' 

For each R > consider the compact set i^^ := [-R, R] \ Uj ^(zj). By (16.201) . to each 
point x G we can associate an open neighborhood A^(x) that should be sufficiently 
small so it does not contain any of the Zj and 

\u\t,y) -u(t,x)\ + \u(t,y) -u(t,x)\ < 77, y e N v (x) (6.24) 

holds. By compactness we can extract finitely many points z m +i, . . . z p so that the whole 
family of N v (zj) form a covering of [-R, R\: 

p 

l-R,R}c{jN v ( Zj ). 

j'=i 

Note that, in view of the Lipschitz continuity property of the map $ (see (16. 2ft ). the 
graph distance associated with u e and u is clearly controled in the L°° norm. Such an 
estimate is used within each region of small oscillations, as follows. Set uf ,h := u e (t, z|±) 
and := w(t, Zj±) for j = 1, . . . , m. By distinguishing between regions of large jump and 
regions of small oscillations we can estimate the graph distance between the approximate 
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graphs (X e , U e ) and the ^-completion (Y, V$), as follows (where the constants Ci, C 2 , . . . 
depend upon Cip(Q) only) 

dist((X £ , U e ), (Y, V*)) < D 1 + D 2 
with, for the shock regions (J — 1, . . . m), 

£>i = 2 sup sup |w e (i, y) — u(t, x) | 

j=l,...,m x€N v (zj),xf^Zj 

+ 2£ip($) sup sup |w(t, y) — w(t, x)| 

+ 2£?p($) sup sup |w e (t, y) — w e (t, x)| 

j'=l,...,ma;,j/eJV„(aj),x,j/^j 

+ sup sup \§(s]u~ ,h ,up h ) — §(s]Uj ,11^)1 
j'=l,..,m ( e[0,l] 

and, for the regions with small oscillations (j — m + 1, . . .p), 
D 2 = sup sup \u e (t, y) — u(t, x)\ 

j=m+l,...,p x ,y&N v (zj) 

+ £ip(<3>) sup sup \u(t, y) — u(t, x)\ + \u e (t, y) — u e (t, x)\. 

j=m+l,...,p x ,y£N v {zj) 

Using again the Lipschitz continuity of $ we can write 

ms;uj' h ,up h ) - S>(s;uj,ut)\ < Cip{<$>) (\u~' h -uj\ + \u+> h - u+\). 
Therefore in view of (16.231) and (I6.24p we find 

dist((X e ,[/ e ),(Y,^)) < Crj. 
where the constant C only depends upon Cip(<&). This completes the proof of Theorem 16. 11 
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